GLOBAL WELL-POSEDNESS AND STABILITY OF 
ELECTRO-KINETIC FLOWS 

DIETER BOTHE, ANDRE FISCHER AND JuRGEN SAAL 

Abstract. We consider a coupled system of Navier-Stokes and Nernst- 
Planck equations, describing the evolution of the velocity and the con- 
centration fields of dissolved constituents in an electrolyte solution. Mo- 
tivated by recent applications in the field of micro- and nanofiuidics, we 
consider the model in such generality that electrokinetic flows are in- 
cluded. This prohibits employing the assumption of electroneutrality 
of the total solution, which is a common approach in the mathematical 
literature in order to determine the electrical potential. Therefore we 
complement the system of mass and momentum balances with a Poisson 
equation for the electrostatic potential, with the charge density stem- 
ming from the concentrations of the ionic species. For the resulting 
Navier-Stokes-Nernst-Planck-Poisson system we prove the existence of 
unique local strong solutions in bounded domains in R" for any n > 2 as 
well as the existence of unique global strong solutions and exponential 
convergence to uniquely determined steady states in two dimensions. 
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1. Introduction 

Electrokinetic effects have received a lot of attention from various 
branches of research and applications, recently in particular in the field of 
micro- and nanofluidics. Two well-known phenomena are electrophoresis and 
electro- osmotic flow (EOF). While in electrophoresis, a dispersed (colloidal) 
particle is dragged through a fluid under an electrical field, EOF describes 
the motion of an aqueous solution past a solid wall as a response to an ex- 
ternal electrical field. The latter offers a powerful method to manipulate 
liquids in mini- and micro-devices, since it does not depend on mechanical 
tools. Indeed, employing electrical fields to pump a liquid through a (sub- 
) micro-channel can economize this task considerably, if there are charged 
solutes in dispersion; see, e.g., [281 [33]. EOF can be applied, for instance, 
for cooling systems in micro-electronics. In analytical chemistry, electrical 
fields serve for the separation of chemicals according to their electrophoretic 
characteristics (see |24] ) . but it also has a great potential for promoting 
mixing on micro-scales [8]. For further applications we refer to [51] and the 
references therein. 

The present paper studies the mathematical model governing electro- 
osmotic flow. From the significance of the EOF, a proper mathematical un- 
derstanding is indispensable for applications. We consider a rather general 
situation in which a dilute viscous solution with dissolved charged species 
is placed in a container with solid walls, the container being situated in an 
electrical field. A typical concrete case is an aqueous electrolyte solution, 
where of course uncharged species may also be dissolved in the solute. The 
externally applied electrical field induces the EOF due to the presence of 
charged constituents in solution. On the other hand, the ions themselves 
generate an intrinsic electrical field as well, which in turn also affects their 
fluxes. Thus, even without external electrical fields, electrical effects have 
to be accounted for whenever charged species are involved. Concerning the 
solid-liquid interface, it is well-known from electro-chemistry, see e.g. [31j . 
that electro-chemical double layers appear on the surface, thus the surface is 
typically charged. In a simplified view the surface can be interpreted locally 
as a plate capacitor. In contrast, far away from walls or other interfaces, 
the net charge density is essentially zero, i.e. the solution is electro-neutral 
inside the bulk; see, e.g., [221 131] . But note that in nanofluidic applications 
as well as in transport processes in nanopores within porous media, there is 
no macroscopic bulk phase and, hence, electroneutrality does not even hold 
in a subdomain. Even the more general approach to describe the electric 
potential via the Poisson-Boltzmann equation may not be accurate enough, 
since the Boltzmann statistics does not apply in overlapping double layers; 
cf. [30]. 
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In the mathematical literature, the electroneutrality condition has been 
frequently used as a simplifying condition, see e.g. [31 El 02]. As a conse- 
quence of electroneutrality, there also is no electrical body-force acting on 
the mixture, such that the equations for the fluid motion decouple. There 
still is a non-trivial electrical field, implicitly determined by the algebraic 
condition of electro-neutrality which leads to a strong coupling of the mass 
fluxes, but still this simplification is too restrictive here, since it rules out 
the electrokinetic effects. 

In papers of Y.S. Choi and R. Lui, see j9[ [TO], as well as in papers of 
H. Gajewski, A. Glitzky, K. Groger and R. Hiinlich, e.g. [TOJ [TOJ EQl EE] , 
the electro-neutrality condition is not employed. Instead, in the context of 
reaction-diffusion systems, a more general mathematical model comprising 
species mass balances with fluxes according to the Nernst-Planck equation 
combined with a Poisson equation for the electrical potential is investigated. 
However, no fluid motion is taken into account and, hence, no momentum 
balance is accounted for. Employing suitable Lyapunov functionals, the exis- 
tence of global solutions in two dimensions is proven and a careful analysis of 
the long-time behavior is provided. In this approach the mass flux ji of con- 
stituent i is composed of a Fickian diffusion term and an electro-migration 
term due to the electrical field. 

In [35], the afore mentioned system of Nernst-Planck equation and the 
Poisson equation is complemented by the Navier-Stokes system. Well- 
posedness issues are accounted for under the assumption that all diffusivities 
of the individual species are constant and equal. From the physical point 
of view, these assumptions are rather strong. In fact, the diffusivity of one 
species will in general depend on the full composition of the system; cf. 
[TT] . While this dependence is negligible in the dilute case considered here, 
the values of the diffusivities of individual constituents will still differ sig- 
nificantly, up to orders of magnitude. Moreover, in |35] the electrostatic 
potential is assumed to satisfy the homogeneous Neumann condition, which 
rules out boundary charges in contrast to the possible occurrence of double 
layers at the boundary. 

Further in some way related mathematical literature is given by [5[ [23] . 
In [5], for instance, motion of a rigid charged macro-molecule immersed in 
an incompressible ionized fluid is considered and local existence of a weak 
solution is proved. The work [23] deals with a two carriers (anions and 
cations) electrophoretic model in the whole space M. n . There local-in-time 
existence of a unique strong solution in a Hilbert space setting is proved. A 
major point in [23] is uniformity of the existence interval with respect to the 
viscosity, allowing for the inviscid limit to the Euler-Nernst-Planck-Poisson 
system. 

Apparently, no further analytical works are available concerning the sys- 
tem consisting of the Navier-Stokes equations for the fluid motion, the 
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Nernst-Planck equation for the species concentrations and the Poisson equa- 
tion for the electrostatic potential. Therefore, our approach generalizes the 
ones given, for instance, in [9l QUI CG3 HU EUJ [21] , where electro (-reaction)- 
diffusion systems are considered, not taking into account convection induced 
by fluid motion or a coupling to the Navier-Stokes equations. In this con- 
text, we also want to stress the fact that our approach does not need the 
electro-neutrality assumption, as is assumed e.g. in [3JE1I32], and also allows 
for different diffusivities for the individual species (in contrast to [35]) as well 
as for non-trivial charge density on the interface. Our result on local-in-time 
well-posedness for the full Navier-Stokes-Nernst-Planck-Poisson system in a 
strong L p -setting hence is new. The methods used for the proof of global 
strong solutions in two space dimensions are close to the ones utilized in 
[10] , The results obtained in the underlying note, however, differ from the 
ones derived in |10] in several respects. Apart from the fact that in [10] 
there is no convection, an approach in Holder spaces is presented whereas 
we prefer to work in an L p -setting. Our proof of existence and uniqueness of 
steady states is much simpler and we are able to remove a restrictive global 
electro-neutrality condition imposed in [10] for uniqueness of steady states. 
Moreover, we prove exponential stability of steady states which is not given 
in [10]. The methods applied here are related to [20] , where exponential 
stability of steady states for concentrations is derived in the absence of fluid 
motion. In this regard, we note that the proof of exponential convergence to 
a steady state simultaneously for fluid velocity and concentrations of species 
is new. 

This article is organized as follows. First, we give a justification of 
the mathematical model considered in this work, the Navier-Stokes-Nernst- 
Planck-Poisson system, by illustrating its derivation from basic principles. 
We then state our main results on local well-posedness in arbitrary dimen- 
sion and on global well-posedness and exponential stability in two dimen- 
sions. In Section [2] we give a proof of the local-in-time well-posedness result 
based on maximal regularity. Utilizing an appropriate Lyapunov functional 
and conservation of mass for the concentrations, which gives the starting 
point for a bootstrap argument, we prove in Section [3] that in two dimen- 
sions the local solution extends to a global one. The content of Section U] is 
the proof of existence of a unique steady state solution to the Navier-Stokes- 
Nernst-Planck-Poisson system for each given distribution of initial masses. 
Finally, in Section [5] we prove exponential stability of this stationary solu- 
tion. 

1.1. The mathematical model. The basic continuum mechanical model 
for transport processes in (non-reactive) electrolyte solutions consists of the 
Navier-Stokes equations for the mixture as a whole, coupled to a set of 
species equations for which the mass fluxes are modeled via the Nernst- 
Planck equation with the electrical potential given by a Poisson equation; 
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see [3T]. Then, depending on the concrete situation, the Poisson equations is 
often replaced by either the above mentioned condition of electro-neutrality 
or by the Poisson-Boltzmann equation in which the concentrations are re- 
lated to the electrical potential by a Boltzmann statistics. Here, in order 
to maintain the applicability to EOF, we keep the more general Poisson 
equation. 

But there is one hidden difficulty which forces us to briefly reconsider the 
derivation of the mathematical model. Indeed, the Nernst-Planck equation 
in the usual form, stating that the diffusive (molar) mass fluxes ji are given 
as 

F 

ji = -DiVci - mi—ZiCiX7(p, (1.1) 

already relies on additional assumptions. In fact, the standard derivation 
from the more fundamental generalized Maxwell-Stefan equations ([26]) al- 
ready considers the mixture to be electrically neutral. But the resulting 
form is also applied to the general case which obviously requires clarifica- 
tion. Fortunately, it turns out that the constitutive equation (jl.ip is also 
valid without electro-neutrality if the electrolyte solution is dilute - the case 
we are considering here. 

We consider a fluid composed of N + 1 constituents with molar concen- 
trations Ci, i.e. Cj = Pi/ Mi with pi the mass density and Mi the molar mass, 
where i = refers to the solvent in which the other species (i = 1, . . . , N) 
are dissolved. A sufficiently fundamental starting point is the system of 
partial mass balances together with a common mixture momentum balance, 
where the diffusive fluxes of species mass are determined from the system of 
Maxwell-Stefan equations. The latter follow as a reduced set of balances of 
partial momenta, simplified by neglecting diffusive waves due to differential 
acceleration of the constituents. The considered balances read 

dtci + div(au + ji) = 0, i = 0,...,N, (1.2) 
d t (pu) +div(pn® u - S) = -Vvr + pft, (1.3) 

where 

N N N 

p:=^2 Pi' P u := ^2 P iUi " P b := Pibi 
i=0 i=0 i=0 

define the total density p, the barycentric velocity u and the total body force 
b, given the individual velocities and forces 6j. The molar mass flux ji 
is defined as ji := Cj(uj — u), but also has to be modeled, since individual 
velocities are not accounted for. In (|1.3|) . S denotes the (viscous) extra stress 
and 7r the pressure. Note that equations (|1.2|) are not independent, since 
total mass is conserved. Indeed, 

N 

^2Miji=0, 

i=0 
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and one of the species equations in (ll.2p can be replaced by the mixture 
continuity equation 

dtp + dw(pu) = 0. (1.4) 

The closure of (jl.2p and (|1.3j) . i.e. the modeling of the ji is achieved by the 
Maxwell-Stefan equations. Formulated for molar mass concentrations and 
under isothermal conditions they read (cf. [26]) 

N 
fc=0 

Here Xj = Cj/c with c := S^ C/t are the molar fractions, £>ij the Maxwell- 
Stefan diffusivities, R is the universal gas constant, T the absolute tempera- 
ture and yi := pi/ p the mass fractions. For extension to the non-isothermal 
case, see e.g. Chapter 24 in (gj. The molar-based chemical potentials pi are 
defined as 

d(pip) d(fnp) 



Pi = M, 



dpi dci 

where pip = pip(T, Po, • ■ ■ , Pn) denotes the free energy density of the mixture. 
Hence pi are functions of T and all pi as well, but in the literature on 
chemical thermodynamics and multicomponent transport, it is common to 
use (T, it, Xq, . . . , xjv) as independent variables instead. Of course not all 
independent, hence (T,ir, x±, . . . , xjy), say, is more appropriate. If 
Pi = Pi(T,n, Xj) denotes this function, then (jl.5p can be written as 

N ■ - ■ rh - 

- Y Xk3t XiJk = ^V n T pi + ^^Vtt -lLf bi -b), (1.6) 
^ £> ik RT n ' lfJJl ^ RT RT h v ; 

A;=0 

where <j>i is the volume fraction of species i and 

\ - dpi 

j 3 

is the "gradient taken at constant ir and T" . The three terms on the right- 
hand side of (jl.6p correspond to diffusional transport due to composition 
gradients, pressure gradients and external forces, respectively. 

In case of an electrolyte solution, the individual body force densities are 

bi = -^t^V^ (1.7) 
Mi 

where F is the Faraday constant (charge of one mole of electrons), %i are 
the charge numbers and 4> is the electrical potential. The latter is assumed 
to be determined by the electrostatic relation 

N 



div(eV0) =Fj2z k c k , (1. 



k=0 
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i.e. dynamic effects and the influence of magnetic fields are ignored which 
is usually sound for the applications which have been mentioned above. In 
(|1.8p . the material parameter e denotes the electrical permittivity of the 
medium. 

The system of equations ([0]) . (TOD . (TQjl . (fLol) . (TOD and ([0]) has to be 
complemented by an equation of state, relating ir with (T, po, . . . , pjy), and 
by a material function modeling the free energy of the mixture; alternatively, 
and more common, the chemical potentials need to be modeled. These are 
highly complicated tasks, in general. 

From here on, we only consider the case of a dilute solution, for which 
species 0, say, is the solvent with xq ~ 1 and all other species are present 
in trace quantities, i.e. ^ <C 1 for i = 1,...,N. This leads to several 
simplifications. First, for i = 1, . . . , N, the Maxwell-Stefan equations reduce 
to 

Next, the difficult request for an equation of state is resolved by assuming 
the mixture to maintain constant density, which is reasonable since p ~ po 
and the latter will be almost constant at isothermal conditions. Therefore 
the species equation for i = is omitted, but 

p = const, divu = (1-10) 

is used instead. The effect of pressure diffusion is usually negligible, since 
4>i/yi is not far from 1 and here both and j/j are even close to zero. It 
remains to model the chemical potentials. For a dilute solution one has, [3], 

ft = /#(r,7r) + i2riogXi, i = l,...N. (1.11) 

With these simplifications and with Di := DiQ, the diffusive fluxes ji from 
()1.9P become 

3i = -D i (cVx i + ^{z i c i -^z k c k )V^\, i = l,...,N. (1.12) 

This is close to the Nernst-Planck form, up to the term Ylk=o z k°k- The 
latter vanishes in the electro-neutral case, but this would restrict the ap- 
plicability of the model too far. Instead, we note that in the dilute case 
necessarily zq = 0, hence 6o = 0; and then 

N N 

bi-b = bi-^2 Vkh = (1 - Vi)h + ^2 y^ b k i = l,...,N. (1.13) 

k=0 i^k=l 

Consequently, b is small as compared to 6j (for % such that Z{ ^ 0), but 
this does not imply b to be small as compared to the other quantities in the 
momentum balance. In fact the electrical forces acting on dilute components 
can drag the whole mixture as exploited in EOF. To sum up, using the fact 
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that c « Co is nearly constant, we see that in a dilute electrolyte solution the 
diffusive molar mass fluxes of the solutes can be modeled via the Nernst- 
Planck constitutive equation 

ji = -Di {Vci + -^zmvA , i = 1, . . . , N. (1.14) 

Note that this derivation from the Maxwell-Stefan equations automatically 
also yields the Nernst-Einstein relation, saying that the mobility coefficients 
in the electromigration term, i.e. the mi in (jl.ip equal Di/RT. 
Concerning (|1.3p . the viscous stress tensor S is modeled by 

S = fi{Vu + (Vu) T ) 

with constant viscosity fi > 0, hence divS" = fiAu since u is divergence 
free by (jl.lOp . The permittivity e in equation (jl.8p for the electrostatical 
potential is assumed to be positive and constant. Thus the full model reads 

JV 

p{p t u + (u ■ V)u) - fiAu + V7r = —F ZjCjVcj), (1.15) 

3=1 

divn = 0, (1.16) 
d t ci + div (au - Di (va + -^z iCi V<^j \ = 0, i = 1, . . . , N, (1.17) 

AT 

-eAcj) = F^zjCj, (1.18) 
i=i 

for (t,x) G (0,T) x f2, where C K n , n > 2, is the domain under considera- 
tion. We note that the force term on the right hand side of (|1.15p is usually 
called Coulomb force. 

At dO,, we impose the boundary conditions 

u = 0, (1.19) 
-A(&cj - ^ztCiW) = 0, i = 1, . . . , N, (1.20) 

$,0 + T0 = f, (1.21) 

for (t,x) G (0, T) x 50, where ^ denotes the outer normal to O and c^, 
denotes the derivative in direction v. 

Condition (|1.19p is the commonly assumed no-slip condition for the ve- 
locity field u and the no-flux condition (jl.20p models impermeable walls and 
assures the conservation of mass in Q. Relation (|1.2ip requires some physical 
explanation. By the Maxwell equations we have 

ed v cj) = a s on <9f2, (1-22) 

where a s denotes the surface charge density. A typical choice for boundary 
conditions for the electrical potential (see e.g. |10^ I35j) are given by homo- 
geneous Neumann conditions, implicitly claiming that the surface is charge 
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free. As already mentioned, this condition is not sufficient for our situation. 
However, finding appropriate models for the charge density on the surface is 
a delicate matter. The boundary condition (|1.2ip is the same as employed 
in |20j . For a physical motivation we imagine to append a layer of (dimen- 
sionless) thickness < 6 <C 1 around the volume which gives us locally a 
plate capacitor for which the surface charge density is known to be 

cr s = ^{<t>5 ~ 4>), 

where (f>g is the electrical potential at the outer boundary of the layer. Plug- 
ging this equality into (|l.'22p we obtain 

1 1 

d v <j> + -(f> = -(f) S . 

The parameter r := 1/5 > in (jl.2ip hence refers to the local capacity of 
the boundary. The right-hand side £ in (jl.2ip depends also on an external 
electrical potential. 

The system of equations (|1.15p - (|1.18p . (|1.19jl - (jl.2ip complemented by the 
initial conditions 

u(0,x) =u°{x), i£Q, (1.23) 

c(0,x) = c°(x), ieO, (1.24) 

represents the full model, where here and in the sequel c = (ci, . . . , cjv) 
always denotes the full vector of concentrations. 

1.2. Main results. Before introducing our main results let us fix some 
notation which will be frequently used below. 

Let n G N and Q C MJ 1 be an open domain with boundary dQ and 
1 < P < oo. By LP(n),LP(dn),W m >P(n) we denote the usual Lebesgue 
and Sobolev spaces. We do not distinguish between spaces of functions and 
spaces of vector fields, i.e. we write also L p (0) for L p (fi) n for example. The 
notion W™ ,p {Si) describes the closure of smooth and compactly supported 
functions Cg°(fi) in W m ' p {Q). For s > 0, p e (l,oo) we write W°{n) = 
(Lp(£l),W m ' p (£l))Q !P for the Sobolev-Slobodeckii space, where (•, -)g >p is the 
real interpolation functor with exponents 8,p and s = 9m. It is worth 
mentioning that W^{Q) = W m > 2 (n) for m G N, see @Q|. Let C^ a {fl) 
denote the solenoidal functions with compact support in £1 and 

denotes the space of solenoidal L p -functions. Occasionally we write W^^ift) 
instead of W™ (n) n L&(Q). The corresponding norms are denoted by 
II ' ||p: II • \\p,dU, II ' ||vy m -P) etc. For J = (0, T), T e (0, oo] and X a Banach space 
we write L p ( J; X), W m,p ( J; X), W p ( J; X) for the corresponding spaces of 
X- valued functions. For s ^ 1/p let 
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Note that W£(J; X) = W°(J; X) for s <l/p and that for 1/p < s < 1+1/p 
this space coincides with all functions in Wp(J, X) having vanishing time 
trace at zero. 

Now let f2 C W 1 , n > 2 be a bounded smooth domain. Concerning the 
model derived in Subsection 11.11 we set the parameters p, p, F, RT all to 1 
for technical simplicity; all the results obtained in this work remain valid in 
the general case where p, p, F, RT > 0, but constant. 

Remark 1. From the functional analytic approach to the Navier-Stokes 
equations it is sufficient to solve for the velocity field rather than both for 
the velocity and the pressure field, because the pressure can be recovered 
once the velocity is known, see e.g. |37| . In this spirit, by formal application 
of the Helmholtz projection P (see [37]) to (|1.15p . solving system (|1.15p . 
(|1.16p . (|1.19p . (|1.23p is equivalent to solving the problem 

d t u + A s u = -P{u-V)u- P(Y,jZjCjV<t>) , <>0, (1.25) 

u\ t =o = u°, (1.26) 

where ^4^ = —PA is the Stokes operator subject to homogeneous Dirichlet 
boundary conditions with domain 

v{A s ) = w 2 *(n) n w^(n) n i%(n), P e(i,oo). 

Once (jl.25p . (jl.26p is solved, the pressure can be recovered by 

Vtt = (J - P) (Au - (u ■ V)it - (ZjZjCjVcj))) . (1.27) 

Remark 2. Consider problem (fTTT8|) . ([L2T]) . 

(a) For the Robin-Laplacian — A/j in L p (tt) with T>(— An) = {v £ 
W 2,P (Q.) : d u v+rv = 0}, r > a constant, it is well-known that £ p{— Ap>), 
therefore the operator (— Ar,B)v := (—Av,d u v + tv), is a bijection from 
W m + 2 >P(n) to W m > p (tt) x W™ +1 ~ 1/p (dn), m £ N and 1 < p < oo, see e.g. 
@0l Chapter 5]. 

(b) We only consider time-independent functions £ = £(x). This will 
prove important in various places throughout this work. 

(c) Note that once c is known, the potential 4> is known as well. 

From Remarks [1] and [2] it suffices to consider solutions (u, c) rather than 
solutions (u, c, 7r, <j)) of the following problem: 

d t u + A s u + P(u ■ V)ii + P (EjZjCjVc/)) = 0, t>0, (1.28) 

d t Ci + div (au - DiVa - A^CjV0) = 0, t > 0, x £ Q. (1.29) 

-eA0 = XjZjCj, t > 0, x £ ft (1.30) 

c^Cj + ZiCid u (j) = 0, t > 0, x £ <9f2, (1-31) 

d u( f) + T(j) = C, t>0,x£dn (1.32) 

u(0) = u°, x £ n (1.33) 
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cj(o) = c°, xen, (1.34) 

for i = 1, . . . , N. System (fOHjl - ffOIjl will be called problem (P). 
We will frequently use the abbreviation 

Ji ■= ji + CiU = -DiVa - DiZiCiV<j) + CiU. (1.35) 

With this notation and the no-slip condition for u, equations (|1 .29|) and 
(jl.3ip can be written as 

dta + div Ji = 0, in(0,T)xO, i = l,...,N, 
Ji-u = 0, on (0,T) xffi, i = l,...,N. 

For the local well-posedness in any dimension n > 2 we apply maximal 
regularity. Let us define the usual spaces of maximal L p -regularity for 1 < 
p < oo by 

Kr, P ■= W 1 ^(0,T;TJP(n))nlJ'(0,T;V(A s )), 
K c T<p := W X ' P (0, T- L P {£1)) n L p (0, T; W 2 ' p (n)), 
E TtP := E^ p x E^, 

x p ■.= (w^nn) n w$*(n)) x ^ p 2 ~ 2 /^). 

We also set 

:= L p {<d,T-W 1 ' p {0)) 

E^ p := W^O, T; W 2 ' p (n)) n L p (0, T; ^ 4 ' p (0)) 

for the spaces of the corresponding pressure and potential. Our main result 
on local-in-time strong solvability now reads as 

Theorem 1.1. Let £1 C W 1 , n > 2, be a bounded and smooth domain and 
let (n + 2)/3 < p < oo, p / 3. Let £ € Wp~ 1/p {dn) and {u°,c ) e X p . In 
the case p > 3 let c° and £ satisfy the compatibility condition 

d„$ + z4d„<t>°\ dU = 0, (1.36) 

w/iere 

T/ien i/iere is T > such that there exists a unique solution (u, c) E Et, p to 
problem (P). The latter is equivalent to saying that 

(u,c,7r,0) eEr iP x^xE^ 

wii/i 7T given by {1.27\ ) is the unique solution of system U.15\) - [1.21\) and 
fl.23\) . {1.21$ . If c° > (componentwise) we have c > as long as the 
solution exists. 
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Remark 3. (a) The local-in-time solution (u, c) remains unique, as long as 
the solution exists, see Remark UJ 

(b) Remark [2] and the fact that dt and A commute cause the potential <j) 
to be contained in the following regularity class: 

<f> G W 1 ' 13 ^, T; W 2 > p (fl)) n L p (0, T; W 4 ' p (fl)), 

if £ is smooth enough. This is accomplished by £ G Wp ^^(dQ) which 
motivates the assumption on £ in Theorem II. 1[ Note, however, that this is 
done for simplicity; this constraint can be weakened, while still obtaining 
the same results presented in this paper, see Remark [SJ 

(c) We also remark that at this stage minimal boundary regularity is none 
of our major objectives. By the higher regularity for <fi we work with, the 
approach presented in this note should at least work for domains O of class 
C 4 . This is also more than enough to guarantee maximal regularity for the 
Stokes equations, see e.g. [32] , A closer inspection of what is really needed 
for the potential 4> hi the proofs of our main Theorems (which is not done 
here) and again [32] will yield that even C 3 shall be enough. 

(d) The compatibility condition (jl.36p is required for p > 3 only, because 
the trace in time for t = of the boundary spaces is only well-defined in 
this case. 

(e) The main difficulty in proving Theorem II .11 lies in the fact that (jl.3ip 
represents a nonlinear boundary condition which has to be treated in the 
according trace space 

(f) From the physical point of view it is important to note that the total 
mass of each species is conserved due to the no-flux condition at the bound- 
ary; by nonnegativity this implies ||cj(t)||i = ||c?||i as long as the solution 
exists. 

Now, let (u, c) denote the local strong solution for the case p = 2 and let 
the total free energy function E be defined by 

1 N f e et 

E(u,c) :=^\\u\\ 2 2 + J2 / c i logCidx+|||V^||| + y|M||, 9n; 

for (u,c) G ^(Q) x L 2 (0) with c; > 0. Here ^||n||2 represents the kinetic 
energy, f Cjlogq refers to the chemical potential for species i, see e.g. [HI] , 
and the last two terms correspond to electrical energy. It turns out that the 
function V, given by 



V(t) := E(u(t),c(t)), t>0, 
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is a Lyapunov functional for system (P) and its derivative 4^V, the dissipa- 
tion rate, is given explicitly by 

j N 

-V(t) = -\\Vu(t)f 2 - / jr^simfdx < 0, (1.37) 

where ji = —D^SJci — DiZiC{V(j) is the flux of species i. For a suitable 
definition of the possibly degenerate integral we refer to Lemma 13.71 

In the two-dimensional case this is the cornerstone for the following result 
on global well-posedness. 

Theorem 1.2. Let C K 2 be bounded and smooth and let (u°,c°) € 
Wq£(Q) x (W 1 ' 2 ^) nL°°(0)) with c° > 0, i = 1,...,N. Then the lo- 
cal strong solution (u,c) £ E^2 to (P) (equivalently, (u, c, 7r, (j>) G Et,2 x 
E t,2 x E t,2 t0 (LW-mW> P^W ) from Theorem^ extends to a 

unique global strong solution. 

Remark 4. (a) Comparing the requirements on the initial data of Theo- 
rem ll.ll for p = 2 and of Theorem ll.2l taking into account W^ity = W 1,2 (Q) 
we note that there is an additional L°°-condition on c . This is inserted for 
technical reasons. However, it is possible to ease this constraint, as explained 
in Remark 1131 

(b) We emphasize that for the global existence of solutions due to Theo- 
rem 11.21 no smallness or compatibility conditions on the data are required. 
We also would like to stress at this stage that concerning 3D not only the 
Navier-Stokes but also the Nernst-Planck equation seems to be difficult to 
handle. So, even under the assumption of small data, to the authors at 
the present stage it is unclear how to extend Theorem 11.21 to three space 
dimensions, see also Remark 1111 

(c) From the fact that V is a Lyapunov functional we readily confirm 
that the quantities ||ti||2 arid 1 1 log Cj 1 1 1 are uniformly bounded in time. 
Having this information at hand it is possible to prove that the L 2 (f2)-norm 
of the concentrations Cj is uniformly bounded in time, as it is done in [10j . 
Applying further energy estimates, the norm of the strong solutions (u, c) is 
shown to remain finite for finite time, whence the solution can be extended 
up to any time T. 

In view of the long-time behavior of system (P) we are interested in 
steady state solutions in two dimensions, i.e., in existence and uniqueness of 
solutions (u 00 ^ 00 ) to the stationary system 

A s u + P{{u ■ V)u) + P(EjZjCjV(t>) = 0, 

div(cjii — DiVci — DiZiCiVcp) =0, x G Q, i = 1, . . . , N, 

—eAcf) = T,jZjCj, x £ Q, 

dyCi + ZiCid u (p =0, x G 30, i = 1, . . . , N, 

d v 4> + T 4> = f > x ^ dfl, 

(Ps 
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subject to cf ) {x)dx = rrii for given total masses rrii > 0. 

Theorem 1.3. Let $7 C M? be bounded and smooth and let numbers rrii > 0, 
i = 1, . . . , N, be given. 

(a) Then there is a unique solution 

KV 00 ) g (wgCfi) n w 2 ' 2 (n)) x w 2 ' 2 (n) 

to problem (P s ) subject to cf° > 0, J^cfdx = m%, i = l,...,iV. More 
precisely, we have u°° = and c°° G W 2 ' 2 (tt) is the unique nonnegative 
solution to 

div(- AVcj - DiZiCiVj) = 0, i = !,...,#, 



(c^c* + ZjQ^^lan =0, i = 1, . . . ,N, 

—eA<j) = J^jZjCj, 
(d v (j) + T(f>)\dn = £, 



(Pi) 



subject to f n cfdx = m^, i = 1, . . . , N. 

The pressure 7r°° in the equilibrium state is given by 7r°° = Ej-cS up to a 
constant. 

(b) Let (u°,c°) G Wffltl) x (W 1 - 2 ^) n L°°(0)) with c° > and 
f n c®(x)dx = rrii, i = l,...,iV. Lei c°° G VF 2 ' 2 (r2) be the unique non- 
negative solution to (P' s ) subject to f n cf°dx = rrii, i = 1,...,N, and let 
(u, c) be the unique global strong solution to (P) with initial data (u°,c°) 
from Theorem \1.2\ then 

lim(n(t),c(t)) = (0,c°°) (1.38) 

t— >oo 

in W^ 2 {fl) x w^ 2 (n). 

It is worth to mention that there are no equilibrium states with non-trivial 
velocity, i.e. with Theorem II .31 the long-time dynamics of the mixture is de- 
termined by the Nernst-Planck equation and the Poisson equation, problem 



( P' s ) , and the initial masses rrii = j c° i {x)dx > 0. 

The proof of Theorem 11.31 can be outlined as follows. For existence we 
employ the global strong solution from Theorem 11.21 Indeed, let < c? G 
W x ' 2 (n) n L°°(Q) with / n c° i {x)dx = rrii, i = 1,...,N, and u° G W^ 2 (Q) 
be arbitrary and let (u, c) be the corresponding global strong solution from 
Theorem 11.21 The estimates derived in the proof of Theorem 11.21 below 
imply compactness of the semi-orbits (it, c) in L 2 (£l) x L 2 (Q), where any 
accumulation point turns out to be a steady state. From mass conservation 
we know that each constituent's mass is conserved, hence a steady state for 
masses rrii is found. 

With the help of an alternative representation of the dissipation rate 
(jl.37p it is possible to obtain a characterization of steady states, which 
allows us to conclude the uniqueness of steady states subject to prescribed 
masses m,- > 0. 
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With this reasoning, assertion (6) of Theorem 11.31 essentially follows as a 
by-product of (a): by compactness of semi-orbits and the uniqueness of the 
accumulation point, the steady state, the system converges to an equilibrium 
as time tends to infinity. This is first obtained with respect to the L 2 - 
topology, but with additional effort convergence with respect to W 1,2 can 
be proven. 

Concerning the rate of convergence, we have the following result. 

Theorem 1.4. Let VL C M 2 be bounded and smooth, let (u°,c°) G Wq'^(Q) x 
(W 1 ' 2 ^) n L°°(fi)) with 4 > and j n c®dx = m i; i = 1, . . . , N, let c°° 6 
W 2,2 {yi) be the unique nonnegative solution to with J^c™ = m^, i = 
1, . . . , N, and let (u, c) be the global strong solution to (P) from, Theorem \1.2\ 
with initial data (u°,c°). Then there are constants C, uj > such that 

||«(t)||2 + ||c(t)-c~||i<C'e- wt . 

In order to derive exponential convergence results, the so-called entropy 
method has been widely used, c.f. |13j . [20] and references therein. The idea 
is the following: Consider a Lyapunov functional ^(ip(t)) with 

where ip represents the state of the (physical) system, the dissipation rate 
D is nonnegative and *$>('ifj) = if and only if ip is an equilibrium. Suppose 
that ^(ip(t)) < CD(ip(t)) holds for some constant C > 0. Then Gron- 
wall's inequality implies \& (?/>(£)) < \£(0) exp(— i/C), i.e. the value of the 
Lyapunov functional decreases with exponential speed. In mathematical 
works it seems quite common to refer to VP as the entropy. Note, however, 
that this notion can be misleading, since the entropy is a physical quantity 
in thermodynamics which is a priori not connected with a generic Lyapunov 
functional. 

In (|5.ip we define 

*(t) :=E(u(t),c(t))-E(0,J>°). 

From the physical point of view vE'(t) can be interpreted as the difference 
of the total free energy at time t and the total free energy in the equilib- 
rium state. The main task for proving Theorem 11.41 lies in the proof of the 
exponential decay of M r . To this end, note that Aty = gjV, so there is an 
explicit representation for the derivative. It is therefore sufficient to show 
*f?(t) < C4r i f>(t) for some constant C > 0. Performing a similar approach 
as in [20j this can be shown by an indirect (somewhat technical) reasoning, 
see Theorem 15.11 The assertion of Theorem 11.41 then follows by a straight 
forward calculation. 
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2. Local in time well-posedness 

Let ft C K n be bounded with smooth boundary, n > 2. By means 
of maximal L p -regularity we will show local well-posedness simultaneously 
for system (OKI) - (Oil and for system (P5]) - (jTM]) with (fL29j) and (Ol) 
replaced by 

d t Ci + div (-AVcj - DiZicfVj) + qu) =0, t > 0, x G ([09T ) 

S^Cj + Zicfd u (f> = 0, t > 0, x G dft, (PIT ) 

where / + (x) = max{0, /(a;)} for / : M™ -> R. 

This auxiliary problem will be denoted by (P')- ^s well-posedness will 
give us the possibility to derive nonnegativity of concentrations Cj in an 
easy way without applying the strong maximum principle. In the case of 
classical solutions (tt, c) the maximum principle implies that concentrations 
Ci are nonnegative, see e.g. [10]. The reason why we hesitate to follow this 
line lies in the fact that, in general, strong solutions (u, c), see Definiton l2.1[ 
are not classical. Nevertheless it is possible to show that (u, c) is classical 
via maximal regularity and bootstrapping. Since due to the strong coupling 
of the differential equations of (P) this is rather technical, we believe that 
our approach to nonnegativity, inspired by [H], is appropriate. 

Definition 2.1 (Strong solution). Let ft C W 1 , n > 2, be bounded and 
smooth, T > and 1 < p < oo. The function 

(u, c) G W lj? (0, T; LP (fi) x L P (Q)) n L p (0, T; V(A S ) x W 2 ' P (Q)) 

is called strong solution to (P) resp. (P') if (u, c) satisfies (P) resp. (P') 
almost everywhere. 

Remark 5. In order to show local well-posedness for (P) it is sufficient to 
show that (P') is well-posed. The reason lies in the fact that the (nonlinear) 
terms in (P') have less regularity as compared to those in (P). To be 
precise, for c G W 1,p (0, T; L P (Q)) n L p (0, T; W 2 ' P (Q)) we only know that 
c+ G W 1 ' p (0,T;L p (O)) nL p (0,T;VF 1 - p (O)), cf. [TBI - Hence, once we have 
proved local well-posedness for (P')> the same line of arguments yields local 
well-posedness for the original problem (P). 

Following this approach we first have to examine the associated linear 
problem and prove maximal regularity for the corresponding solution oper- 
ator. In the next step we write the full problem as a fixed-point equation 
and apply the contraction mapping theorem to deduce a unique fixed point, 
the local- in-time strong solution. 

In order to sort the terms in (P') into linear and nonlinear terms, we split 
the potential (p into (pi and fa satisfying 

-eAfa =0 in ft, d v fa + rfa =£ on 5ft, (2.1) 

—eAfa=T,jZjCj in ft, d u fa + rfa =0 on 3ft. (2-2) 
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In view of Remark [2] the functions <j)\ and <j>2 are well-defined, if the cor- 
responding right-hand sides in (12. ip and (12. 2p belong to suitable regularity 



classes. This yields the following nonlinear problem, 

d t u + A s u + P [EjZjCjVfa] = F(u, c), (2.3) 
d t Ci - DiAci = Gi(u, c), i = l,...,N, (2.4) 

d uCi = H t (c), i = l,...,N. (2.5) 

u\ t =o = u°, (2.6) 
Ci\ t =o = cl i = l,...,N, (2.7) 

where the nonlinear terms are defined by 

F(u, c) = -P((u ■ V)«) - P (XjZjCjVfoic)) , (2.8) 
Gi(u, c) = D iZi [Vc+ • (V0i + VMc)) ~ c+A0 2 (c)] - Vc, • tt, (2.9) 

i = l,...,N, (2.10) 

tf;(c) = -Zj [^(^1 + 2 (c))] , i = l,...,N. (2.11) 

The linear problem associated to (|2.3p - (|2.7p reads as 

d t u + A s u + P [EjZjCjVc/}!) = f, (2.12) 

d t ci - DiAa = g h i = 1, . . . , N, (2.13) 

d uCi = hi, i = l,...,N, (2.14) 

u\ t=0 = u , (2.15) 

Ci\ t =o = c°, i = l,...,N, (2.16) 
with right hand sides in suitable function spaces. 



Strong L p -solvability of this linear problem is characterized in the follow- 
ing lemma. 

Lemma 2.2. Let SI C W 1 , n > 2, be bounded and smooth, (n + 2)/3 < p < 
oo, < T < oo and £ £ Wp 1 ^ p (dQ,). Then there is a unique solution 

u E ^(0,T;^(O)) n LP(0,r;2?(A 5 )) =: E^, 

c E Ty 1,p (0, T; L p (fi)) n L p (0, T; W 2 ' p (0)) =: E^, 

to problem $2.12\) - [2.16\) iff the data is contained in the following regularity 



classes 

feLV(0,T;Ll(n))=:W% tP , (2.17) 

geW(0,T;W(n))=:F^ p , (2.18) 

h E W^ 2 -^ 2p (0, T; L p (dn)) n L p (0, T; W}-^ p {dQ)) =: F% 9 p n , (2.19) 

u° E iy p 2 " 2/p (0) n WgCn) =: X£, (2.20) 

c° E W p 2 " 2 / p (0) =: X£ (2.21) 
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and the initial data satisfies the following compatibility condition 

d v c\ = Johi, i = 1, . . . , N, in the case p > 3, (2.22) 
where 70 is the trace operator in time for t = 0. 

Proof. The existence of a unique solution c G to the parabolic problem 
(|2TT3j) . (gum and (I2TTHD iff the data fulfills flZJgp , (jOty . (I23TD and (jZZZD 
is well-known and contained e.g. in [12] . From £ G VF P l ^ p (dVL) we obtain, 
cf. Remark El 01 G Ty 4 ' p (ft), hence V^i G W 3 ' p (0). Note that X is time- 
independent. For equation (|2.12p we compute 

\\(H^(I>i\\lp(LP) < \\ c i\\LP(LP)\\^4>l\\BUC(n)' (2.23) 

due to W 3 ' p (9,) ^ BUC(£l) from Sobolev's embedding theorem, [I]. So 
/ := / — P(SjZjCjV(/>i) G L p (0, T; L^(^)). By maximal regularity for the 
Stokes operator, see e.g. j38} [17) 152] the remaining claims follow. □ 

Let E T , P := E£. p x E^ p , F T , P := F^ p x F^ p x and X p := X% x 

X p . From Lemma 12.21 for every given T > there is a continuous solution 
operator 

L: {(f,g,h,u°,c°) G F TjP x X p : (12T22D is valid} E T , P , (2.24) 

for problem (j2TT2l - (j2TT6|) . Observe that for h G F^ the trace at t = is 
well-defined if p > 3. For we set 

Fg n := ( Wy 2 - 1 / 2p (0, T; U } {dtt)) n L p (0, T; Wj _1 / p (0J2)) 

and denote by oFt, p the space 

F T , p = F^ p x ¥ T>p x Fg n . 

Lemma 2.3. Let n > 2, (n + 2)/3 < p < oo, p / 3 and < I" < T < oo. 
T/ien i/iere is C > independent of T' < Tsuch that 

H L |oF T ,, p x{0} 2 ll£(F T , !p x{0} 2 ,E T , iP ) < C, 
where L is the solution operator from \2.2J$ . 

Proof. We abbreviate Lt ■= L\ Q ^ Tp . From \29\ Lemma 2.5] for any Ba- 
nach space X there is a bounded extension operator £ T : oWp(0,T'; X) — > 
oWp(0, T; X), whose norm is independent of T' < T for all s G [0, 2]. So we 
may extend (f,g,h) G Ft", p to some (f,g,h) G Ft, p with (/, g,h)\{o.T') = 
(f,g,h). Thanks to the fact that Lx(f, g, h, 0, 0)|( 0j t') — L(f,g,h,0,0) we 
deduce 

\\L(f,g,h,0,0)h T ,, p < \\L T (f,g,h,0,0)\W TtP < \\L T \\\\(f,g,h)\\ o¥T ^ 
<C\\L T \\\\(f,g,h)\\ a¥T , p 

for (/, g, h) G oFt', p? where C > does not depend on T' < T by virtue of 
[291 Lemma 2.5]. □ 
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Knowing the mapping properties of the linear solution operator we can 
now treat the nonlinear problem fj2.3[> - (|2.T[) involving the nonlinear terms 
F,G,H, given in (jZBD - flZJT]) . We let N(u,c) = {F(u,c),G{u,c), H(c)) in 
the following and write fa = fa{c) for the solution to (|2.2j) to indicate the 
(linear) c-dependence of fa. The fact that F,G,H map into the "right" 
spaces when p is chosen appropriately is the content of 

Lemma 2.4. Let n > 2, (n + 2)/3 < p < oo and T G (0, oo). T/ien 
iV(u, c) G F TjP /or any (it, c) G E£ x . 

For the proof we have to estimate the individual terms of N(u,c). How- 
ever, these computations are contained in the proof of Lemma 12.61 and are 
therefore omitted here. 

Let (ti°,c°) G X p be fixed from now on. We will prove the existence 
of a unique local-in-time strong solution to (|2.3p - (|2.7p with (n(0),c(0)) = 
(u , c ). This is equivalent to the existence of a unique fixed point of the 
equation 

{u,c) = L(N(u,c),u°,c°). (2.25) 

It will be convenient to split the solution in a part with zero time trace at 
t = plus a remaining part taking care of the non-zero traces. If c G is 
a solution of the Nernst-Planck system, in the case p > 3 due to Lemma 12.41 
by taking trace of H{c) we obtain 

H(c)\ t=0 = H(c°) G W l -*/P(dVL). 

We set h* := e iAsn i^(c°), where A^q denotes the Laplace-Beltrami operator 
on the smooth manifold dVt. By maximal regularity for Aqq, see e.g. [MIES], 
then we have h* G ¥ c ^ and h*(0) = H(c°). 

If p < 3 we just set h* = 0. According to Lemma 12.21 we may define a 
function (u* , c* ) by 

(u*,c*) := L(0,0,h*,u°,c°). (2.26) 
Subtracting (|2.26j) from (|2.25p yields the new fixed point problem 

(«,c) = L(JV'(«,c) l 0,0), (2.27) 

where 

N(u,c):=(F(u,c),G(u,c),H(c)) 

:= (F(u + u*,c + c*),G(u + u*,c + c*),H(c + c*)-h*), [ ' } 

for (u,c) G Et, p := {(«,c) G E T , p : («(0),c(0)) = 0}. 

The advantage of applying the fixed point argument to the zero trace 
space qEt lies in the fact that the embedding constants for embeddings in 
time do not depend on T. This enables us to choose T as small as we please 
without having the embedding constants blow up. 
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Lemma 2.5. Let p 6 (l,oo), I G N, Sj G [0,2], and 6e Banach spaces 
for i = 1, . . . ,1. Suppose there is q G (1, oo) and a Banach space Y such 
that there is a continuous embedding 

l 

f)W?(p,T ;Xi) ^L"(0,T ;y). (2.29) 
i=i 

Then for each T' < Tq there is a continuous embedding 
I 

f) W?(0,T';Xi) ^L"(0,T';Y), (2.30) 
whose embedding constant does not depend on T' <Tq. 

Proof. For T' G (0, To] let Sj,, be the extension operator 

W;(0,T';X) ^ W;(0,T ;X), 

s G [0,2], from \29\ Lemma 2.5], whose norm does not depend on T'; note 
that Sj, in \29\ Lemma 2.5] does not depend on s G [0,2], that is, for each 
s G [0, 2] we have the same extension operator. Let 1Zt' be the corresponding 
restriction operator 

K T > : W;(0, T ; Y) -> W;(0, T>- Y), 
Then the claim follows, since the following diagram commutes: 

f]LioW°>(0,T';X t ) % nf=ioW£<(0,T o ;X) 

Li(0,T';Y) L?(0,T ;y), 
T < T . □ 

Remark 6. Note that, due to the embedding (|2.29p . the embedding con- 
stant for (|2.30p can depend on To. 

For the proof of Theorem 11.11 it will be necessary to show a contraction 
estimate for the map L o N, cf. (|2.27p . This is somewhat technical due to 
the intricate structure of the trace space F^^ 2 . Therefore we prove it as a 
separate lemma. 
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Lemma 2.6. Let n > 2, (n + 2)/3 < p < oo, < T' < T < oo and R > 0. 
Then there are a, C > 0, independent ofT' and R such that 

\\N(u\c l ) - iV(n 2 , C - 2 )|| FT , p < C\\{u\c l ) - (n 2 ,c 2 )|| ET , p X 

<.R + ||u*||l 3 p(0,T';L 3 p) + l|Vn*|| L 3p/2( 0>T /. L 3p/2) 

+ II c *IIl 3 p(0,T';L 3 p) + II c *IIl 3 p/ 2 (0,T';W /1 > 3 p/ 2 ) 
+ ||V02(C*)|| L 3P(O,T';Z, 3 J') + l|V02(c*)|| L 3p/2 (O)T /. w l,3p/2) 

i (2.31) 

+ \\dtC*\\LP(0,T';LP) + || 9t (V<^2 (c*)) ||iP(0,T' jW 1 *) 



+ r to |||c*|| B[/c;(0jT/;H/ 2- 2/P) + ||V(^i|| W i,3 P 

+ ||V02(c*)||sf7C(O,T';BC/C) 

/or (u 1 , c 1 ), (u 2 , c 2 ) G B ^ T , (0, i?), i/ie closed ball centered at in o^T',p 
with radius R. 

Proof. Before estimating the individual terms in N(u l , c 1 ) — N(u~ 2 , c 2 ) recall 
the following estimate which is obtained as a consequence of the mixed 
derivative theorem; for a, (3 G N, p G (l,oo) it holds for s G (0, 1) that 

W a > p (0,T;LP(ft))nL p (0,T;V^' p (ft)) M- W^O.TjW^-'^fi)). (2.32) 

A proof for the case T = oo and the whole space can be found in [36j . 
Applying suitable extensions relation (|2.32j) follows for our situation, see also 
[34j . With (|2.32p and Sobolev's embedding theorem it is straight forward to 
see that 

W 1 ' p {0,T;L p )nL p (0,T;W 2 ' p ) ^ W* 1 (0, T; W 2 ~ 2si ) ^ L 3p (0,T; L 3p ), 

(2.33) 

V(W l > p (0,T;L p ) n L p (0,T;W 2 > p )) Wp 2 (0, T; M^ _2s2 ) 

^L 3 ^ 2 (0,T;L 3 p/ 2 ) (2 ' 34) 

for si = 2/(3p) and s 2 = l/(3p). Of course (i2T32l - (l2T3i]l are also valid for 
T substituted by T'. 

These facts represent crucial ingredients for the proof. Let us, for the 
sake of readability, introduce the following abbreviations: 

^ _1_2 ~ -1-2 » -1-2 

u := u — u , c := c — c , q := c\ — c\ 

and accordingly c + and c^~. The individual terms in N{v},c l ) — N(u 2 ,c 2 ) 
can therefore be written in the following form: 

F(u x , c 1 ) - F(u 2 , c 2 ) = + u 2 ) • Vu + u • V(u* + u 1 ) 

+ E ^« c i " c')V0 2 (c) + CjVMc* + ^))) , (2 ' 35) 
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Gi(u\c l ) - Q(u 2 ,c 2 ) = A^(Vc+(V0i + V0 2 (c* +c 1 )) 

+ V(c* + cf + )VMc) + (c* + c 2 ' + )A0 2 (c) (2.36) 

+ cfAfcic* + c 1 )) + (u* + u 2 )Vdi + u ■ V(c* + cj) 

Hi{c l ) - Hi(c 2 ) = - Zi ((c* + c 2 ' + )cU 2 (c) + c+d„(<Pi + 02(c* + c 1 ))). 

(2.37) 

Concerning the individual summands in F and G we observe that the 
term resulting from the convection term (« • V)ti has the weakest regularity 
properties as compared to the remaining terms in F and G. Since F, G have 
both to be estimated in L p (0,T'; L p ) it is sufficient to consider only the first 
two summands in (|2.35p : the other terms left can be treated in exactly the 
same way. By the continuity of P and Holder's inequality, there is a constant 
C > not depending on T' such that 

\\P((U* + U 2 )Vu + UV( U * + U 1 )) \\ L P(0,T';LP) 

< C((ll n *llL 3 P(0,T';L3p) + II'" 2 ||l 3 p(0,T';L3p))II V ^IIl3p/2( 0jT /. £ ,3p/2) 

+ (l|Vn*|| i 3j,/2(o i T';Z 3 P/ 2 ) + II^ 7 ^ 1 |Il3p/ 2 (0,T';L 3 p/ 2 ))II^IIl 3 p(0,T';L 3 p))- 

(2.38) 

From the fact that u G o&p' p ano - relations (12.33p . (12.34p Lemma 12 . 5 1 assures 
the existence of a constant C > 0, not depending on T', such that 

INIl 3 p(0,T';L 3 p) + II^IIl 3 p/ 2 (0,T';L 3 p/ 2 ) ^ ^ INIe^,, • (2.39) 

So combining f|2.38|) and (|2.39jl we deduce 

||P((u* + U 2 )Vu + £V(u* + U 1 )) \\ L P(o,T';LP) 

< C\\u\\e t , iP (R+ ||^*||l 3 p(0,T';L 3 p) + II Vu *IIl 3 p/ 2 (0,T';L 3 p/ 2 )) 

For further estimation note that, Remark [2l 

||^2(c)||H/l>P(0,T';W 2 >P)nLP(0,T';W 4 >P) < C||c||k=, , (2.40) 

for c G q&ti ti 5 with some constant C > independent of T'. So in performing 
analogous computations for the other summands of F and G, taking into 
account cf < |cj|, we establish 



IK^G)^ 1 ^ 1 ) - (F, G)(u 2 , ^)\\lp{q,t'-,lp) < C\\(u\&) - (u 2 ^)U Tl x 

(R+ ||« |U 3 p(0,T';L 3 p) + IF IIl 3 p/ 2 (0,T';W' 1 > 3 p/2) + ll c * IIl 3 p(0,T';L 3 p) 
+ ll c *llL 3 P/ 2 (0,T';Vy 1 ' 3 P/ 2 ) + II WllL 3 P/2(0,T';I, 3 P/ 2 ) 
+ I|V(/>2(c*)||l 3 p(0,T';L 3 p) + II V<£ 2 (c*) || L 3p/2 ( 0)T / ;W l,3p/2) , 

(2.41) 

for (^^.(^.^G 5^(0,1?). 

We are left with the estimation of (I2.37P in F^f^ . For the treatment of 
this expression let us note that since p > (n + 2)/3 from (I2.32p . Sobolev's 
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embedding theorem, and by choosing 5 < ^(3p — n — 2) such that s = 
1/p + 5 G (0, 1) we have 

G W°(0,T';W*~ 2s (n)) <-> £[/<7([0,T']; ££7(7(0)). (2.42) 

Because of the smoothness of 90 the normal vector field v is BUC 1 (dQ). 
From BUC^dtt) • Z7(<9ft) <^ L p (<90) and BUC l {dD) ■ W l > p {dti) ^ 

W l 'P{dtt) we see, by interpolation, 0Q], BU&idSl) • Wp 1_1/p (dn) 

Wp 1 1 ^ p (dtV); thus we may estimate 

\\Hi{c l ) - Hi(c 2 ) || F c,ao < C||(c* +c-' + )V^2(c)|| LP(0jT; . wl , p)nVK v2 (0ir ,. LP) 
+ C||c+(V^ + V<fe(c* + c 1 ))|| LP(0iT , ;iyl , p)nH/p i/2 (0jT; . LP) , 

(2.43) 

using the fact that F^f^ is the spacial trace space of L p (0,T' ;W 1 ' P ) n 

Wp 1/2 (0,r';LP), see e.g. [EJ. Note that, due to H^) G o^f^, the con- 
stant C in (|2.43p does not depend on T 1 . This can be seen by a similar 
argument as in the proof of Lemma 12.51 

Let us consider the first term on the right hand side of (|2.43p . First we 
estimate the norm in L p (0, T'; W 1,p ). Applying the Holder inequality yields 

=2,- 



< 



i+cS )^Mc)\\lp(0,T';W^p) 

c i llz,3p/2(o,T';W 1 . 3 J>/ a ) + IIl3p/2( 0jT /. H /1,3p/2- ) ) II V(j)2 (c)\\l3p(0,T';W 1 > 3 >>) ■ 



(2.44) 



Next we use the fact that cj(0) = c|(0) = and hence also fa (c)(0) = 0, 
which gives us by applying (12. 33ft . Lemma f2.5l and afterwards relation (|2,40p 
that 

II « + C-' + )V(/>2(c)|| L p (0iT / ;H /i iP) 

< C(\\c*\\ L3p/ 2 {0tT ^ tWl , 3p/2) + R)\\Mc)\k- Tl p (2.45) 

< C(\\C* || L 3p/2( 0)T / ;V yl,3p/2) +ii)||c||E^^. 

1 /2 

Let us now consider the norm in Wp' (0,T';L p (n)). We will not cal- 
culate directly with Sobolev-Slobodeckij-norms. Instead, since <f>2 enjoys 
rather convenient regularity properties, we can afford to use the interpola- 
tion inequality, 



||(c*+^' + )V^ 2 (c)|| VFpl/2(0ir ,. ip) 

< C(T)\\(c* + ^ + )V0 2 (c)||^ 2 (ar , ;LP) ||(c* +c 2 ' + )VMc)Clp(o,T' ]L py 

(2.46) 

The fact that the constant C{T) in (|2.46p does not depend on T' again 
can be seen by applying the extension operator in Lemma 12.51 to (0,T) and 
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the fact that (c* + c 2,+ )Vc/> 2 (c) G W p 1/2 (0,T';L p ). Expanding the norm in 
W X *(Q,T;U>) gives 

ll( c * + C-' + )V0 2 (c)|| H /i,p ( o,T';Lp) = || (<? + C-' + )V^ 2 (c)|| L p (0iT / ;i p ) 

+ \\d t (c* + q' + )V0 2 (c)|| LP(o>T ,. iP) + ||(c* + q' + )a t (V0 2 (c))|| L p (OiT , ;LP) . 

(2.47) 

Note that, since 2 (c)(0) = 0, we have by relations (|2.40|) . f|2.42j) . and 
Lemma 12,51 that 



\\V<h(c)\\BUC(0,T';BUO < C( T )l|V02(c)|| H/s(OT ,.w3-2 S(Q)) 

P (2 48) 

< c"(r)||0 2 (c)|| VK i,p( OiT /. H /2,p) nL p( OiT / ;H /4, P ) < c"(t) ||c||e^, jp - 

This implies for the second term in ([2.47p that 

W d t(c* + C j 2,+ )V(/) 2 (c)|| Z/ p(o,T';Lp) 

< (II^C? ||iP(0,T';iP) + ll^cf ||z,P(0,T';Z,P))l|V^2(c)||Bf/C(0,T';SC/C) (2.49) 

<c(r)(||atc*|| iP(0iT ,. iP) +i?)||c|| E c, ip . 

For the third term in (|2.47p we compute 

||( C * +C^ + )d t (VMc))\\LP(0,T';L P ) 

* _2 - (2.50) 

< (ll c i IIb[/C(0,T';Z,8) + \\ c i llBC/C(0,T';L9))||5 i (V(/) 2 (c))|| Lp(0jT , ;i ,/ :) 



with 1/q + 1/q' = I /p. From [2] it holds 

W lfi (0, T'; L p (fl)) n L p (0, T'\ W 2 ' p (fl)) B£/C(0, T'; W£- 2 /*(n)). (2.51) 

So the term ||c*||s;y<7(o,T / ;L<2) m (12.500 is well-defined if it holds true that 
Wp 2 ^ P (Q) L q . At the same time we have to assure that W l,p {Q) 
L q since c\V0 2 € L P (W 1,P ). The critical cases occur whenp < (n+2)/2. 
In this situation choose q = np/ (n— 2p+2) and q' = npj (2p—2). So, recalling 

p > (n + 2)/3, it holds that Wp 2 " 2/p (0) and W lj, (fi) L q ' thanks 

to Sobolev's embedding theorem, hence 

\M +C 2 ' + )d t (VMc))\\L P (0,T';L P ) 

< C, (H C ill B; 7C(0,T';H/ p 2 - 2/p ) + ll^ 2 !li?[/C'(0,T^^ 2 - 2/p )^l^^ V ^ 2 ^^ll iJ3 ( . T '; Wl ' p ) 

< C ( r )(H C *ll BC /C(0,T';H/ p 2 - 2 ^) 

(2.52) 

In contrast to the L r - norms in time the appearing BU C-norms of c* will in 
general not become small as T' — > 0. Here we need an additional factor T' a 
which will be provided by estimating the first factor in ([2.460 a little more 
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carefully as in (]2.44p . Indeed, we can estimate 

II ( c * +C 2 i' + )VMc)\\LP(0,T';LP) 

< (ll C i IIl3p/2(0,T';L3p/2) + l|Cj' |Il3p/2( 0>T , ;L 3p/2))||V(/>2(c)|| L 3 P ( 0iT / ;L 3p) (2.53) 

< C(T)T' a (\\c*\\ L3p{Q , Tl . L3p) + R)\\c\\ E o Tip 

for a certain a > 0. Note that here and in the sequel the appearing a's at 
first might be different. At the end, however, we just need one a > and 
we can always choose the smallest of finitely many ones. Thus, as for the 
generic constant C, in the sequel a can change from line to line. 

From relations ([2^6]) . ([2^7)1 . (|2^9l) . (|230D . (12321) . (I233]) . and Young's 
inequality we now obtain 

11(4 + c J 2 ' + )V(/> 2 (c)|| VKp i/2 (0)T ,. LP) 
<C7(T)T to (|| C :|| L 3 P(0iT , L 3 P) + J R) 1/2 

• (ll^cllUno.T';^) + ll c *ll BC /c(0,T';VK p 2 - 2/p ) +i2 ) 1/2 |l £ ll E T', P (2 ' 54) 

< c(r)T /a (||c*|| L 3p (0iT / ;L 3 P ) 

+ II^C*llLP(0,T';Lf) + MWbUC^T'-^- 2 ^) + R )^ E t', p - 

Collecting (|2.45p and (|2.54p then gives 
||( C *+c J 2 ' + )V^(c)|| r ,aa 

T',p 

< C(T)(\\C* \\ L 3p/2( 0tT ,. W l,3p/2) + ||c*||l3p( 0iT / ;L 3 P ) (2.55) 

+ \\9tC*\\ L P(0,T';LP) + T ' a \\ C *\\ B UC(0,T';Wt 2/p ) +R )^ E T>, p - 



Now we turn to the second term in (12.43p . This term is basically of the 
same structure as the first term. For its estimation essentially analogous 
calculations and arguments can be carried out. Thus we will be brief in 
detail here. Also note that V</>i has the same regularity in space as V</>2(c*) 
has, but due to its time-independence it can even be treated in an easier 
way. This is why we leave this term out of the following computation. In 
fact, by merely interchanging the roles of the factors we obtain similarly to 
(g31D and (p^5D that 

\\cf(VMc* + C X ))\\lp{0,T>;W^p) 

< C{T){\\V<t> 2 {c*)\\ L ^ T ,^ P) +R)m\w T , ■ 
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1 /2 

Accordingly for the norm in W p (0,T';L p ), by performing similar argu- 
ments as in (|2.46p - (|2.54p we deduce 



< C{T)T' a y\\SJ(t)2{c*) || BUC(0,T>; BUC) + II ^ <p2 (c*) || L 3 P (0,T';X,3p) (2.57) 

+ ||a t (V0 2 (c*))|| 

LPfO.T'jH 71 -?) + ) I|c||e= • 

By the fact that V0i and V</>2(c*) have the same regularity in space, the 
same terms will appear for V0i. Since V<^>i is independent of time then by 
the Sobolev embedding merely a term as T' a ||V<^i||^i,3 P remains. Hence, 
collecting (|2.56p and ()2.57j) gives us 

||q(V0i + V(fo(c* +c 1 ))\\ c ,an 



T',p 



< C(T) (R + T' a [\\VMc*)\\buc(o,t>;BUC) + ||V0i || 



3p 



+ || V0 2 (c*) ||x,3p(0,TV;W 1 . 3 *') + ll^t( V( ? !, 2(c*))||l,P(0 ) T';Ty 1 .P) ) l|c||E T / ip - 

(2.58) 

Finally, we conclude (12311 from relations (JZ33J, 115351) . (12351) and (I23HD . 

□ 

Having Lemma 12.61 at hand, we are in position to prove Theorem ll.il 

Proof, (of Theorem Recall the solution operator L from H2.24H and 
its restriction to qFx' p - Its operator norm is denoted by ||L||. Thanks to 
Lemma 12.61 we can choose first R > and then T 1 > small enough such 
that 

\\Niu\c 1 ) - N(u 2 ,c 2 )\\ ¥t , p < J-IK^c 1 ) - (n 2 ,c 2 )|| ET , p , 

for (u , c 1 ), (u 2 , c 2 ) G B e t , p (0, R). We may assume that T' > was also 
chosen such that 

||#(0,0)|| F „ =\\N(u*,c*)\\vr, < 



-T',p " V" ' " ^I^T'.p — 2||L|| ' 

since all terms in N(u*,c*) consist of time integrals over given functions. 
Now we can apply the contraction mapping theorem. The estimates 

\\LN{u\c l )-LN{u 2 ,<?)U Tl>p < ||L||||iV(u\cV^V)lk,, p 

<i||(u\cV(^V)||E T> 
for (u 1 ^ 1 ),^ 2 ,^) G B oEt ,JO,R) and 

\\LN(u,c)\\ < ||L||(||iV(n,c) - #(0,0)|| + ||iV(0,0)||) 
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2 L " v J " 2 L 



< R 



for (u, c) G 2? e t , (0, i?) assure the existence of a fixed point and hence 
the existence and uniqueness of a local-in-time solution to (P'). Note that 
our computation implies also the unique local-in-time existence of strong 
solutions to the original problem (P), cf. Remark [5j The claimed regularity 
for the pressure then is obtained as an easy consequence of formula (|1.27l) . 

Let us denote the solutions to (P') and (P) by (u a ,c a ) and (it, c) respec- 
tively and let c® > 0, for i = 1, . . . , N. If we can show that for (u a ,c a ) 
we have c a > almost everywhere, then it is also a solution to (P) and by 
uniqueness we see that (u, c) = (u a , c a ). Let c[ = min{cj, 0}. Note that 

7^H c r~ll! = j n dtC i'~ c i'~ dx = f d t$<h~ dx = J d t c i'~ c i dx - 

Prom the fact that dtd%'~ G L p (0,T'; LP) and cf G E^, it is easy to see 
using Holder's inequality and Sobolev's embedding theorem, that the last 
integral is well-defined. Hence formal multiplication of (jl.29fp by and 
integration over Q gives 

^HcHli + X ( " AVC? " A*# + V0 + c?u o )Vc?'-dx = 0. 
Integration by parts, taking into account (|1.31fp . yields 

J4l|crili + A||V^'-||i = DiZi [ c^V^-Vc^'-dx- [ cfu-Vc^dx. 
z at Jn Jn 

(2.59) 

The first integral on the right-hand side is zero since c"' + • Vc"'~ = 0. The 
second one also vanishes due to solenoidality of u, using 

/ c?u-V4'-dx=l [ u-V(cn 2 = 0. 
Jn z Jn 

Here we used the fact that c"Vc" G L l {Q, T'\ L 2 (Q)), to make sure that the 
integral exists. Finally, integrating (|2.59p over (0, t) gives 

||cr(t)||I + 2A /'||Vcf-||I = ||cf'-||l<0 J 
J o 

for almost all t G (0,T). So c" > almost everywhere, whence (u a ,c a ) and 
(u, c) coincide. This proves Theorem 11.11 □ 

Remark 7. Note that the solution (u, c) is unique up to the maximal time 
of existence T max . This can be seen as follows: Let (u,c),(u,c) G Et, p 
be two solutions with initial value (tt°,c°) G X p , and let (0, T') be the 
maximal time interval such that (u, c), (u, c) coincide, T' < T < T max . 
Suppose T could be chosen in such a way that T' < T. From [2] we have 
(u,c),{u,c) G BUC(0,T';X p ), so therefore (u(T'), c(T')) = («(T'),c(T')) 
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and applying Theorem 11.11 with initial value (u(T'), c(T')) £ X p implies that 
(u,c) and (u, c) necessarily coincide on (T',T' + s), say; this contradicts the 
maximality of T 1 . 

Remark 8. The condition £ £ Wp 1 ^ p (d£l) in Theorem 1 1 . 1 1 can be relaxed 
to, e.g., 

£ £ W p 1+2n/(3p) - 1/p (dfi). (2.60) 
It is a straight forward calculation that with this condition the facts that 
V0i £ BUC(Q) and V 2 </>i £ L 3p (f2) remain valid. Having this information 
at hand, we see that estimate (|2.23p and the reasoning to obtain (|2.58p is 
still true. There are no further points where the regularity of <pi is used. 
However, because of the unphysical dependence on the space dimension, we 
chose to state Theorem II. II in the given form. 

Corollary 1. Let (u,c) be the local strong solution to (P) from Theorem 
IJ.il with initial data (it°,c°) where c® > 0, i = 1, . . . , N. Then the initial 
masses mi := f^c ( -(x)dx are conserved, i.e. 

IkWHi = ||c°||i, *>0, i = l,...N. 

Proof. The claim follows from nonnegativity of Ci and 
d_ 
df 



Cidx = \ dtCidx = — I div Jidx = — / J{ ■ vdx = 0, i = 1, . . . N, 

dfl 



(2.61) 

due to (fL29j) and (fPT|) . □ 

3. Global well-posedness in two dimensions 

Before we turn our attention to the estimation of the local solution (it, c) 
from Theorem 11.11 we collect some known technical results which will be 
employed in the proofs of global well-posedness and stability. Below we 
write 

log + x = max{0, log x}, for x > 0. 

Lemma 3.1. |10^ Lemma 3.2] Let v > be defined on O. Then for any 
p > \, e > and < a < 1, there is a constant C = C(p, a, e, O) > such 
that 

||u||p<e||u(Iog + u) 1 - a ||p + C. 

Lemma 3.2. |10t Lemma 3.3] Let m > 1, q > r > 1 and define 

l _ l 

r q 

7=TTXTT 

n m r 

Suppose that < 7 < 1. Then there exists a with 7 < a < 1 and a constant 
C = C(m, n, q, r, O) such that for all v £ W l,m {yi) with v > 

||«(log + «) 1 " a || 9 < C||i;log + i;||i- a (||Vi;|UM|r 7 + Hi). 
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This result is contained in Lemma 3.3 of [10]. Note that the restriction 
r > 1 assumed there is not needed in the proof. 

Lemma 3.3. Let £1 C W 1 be bounded with smooth boundary. Let 1 < q < oo 
if n = 2 and 1 < q < 2n/(n — 2) if n > 3. Then there is a constant 
C = C(n,q,Q) such that 

\M q <c(\\Vv\\%\\v\\{- a + \\v\\ 1 ), 

for any v G W > 2 (£l), where 



a = XT! 

n "i 2 



and < a < 1. 



This is a direct consequence of the Gagliardo-Nirenberg inequality without 
boundary conditions, for a proof see e.g. [27] , 

For convenience we also state the Poincare inequality in a form which will 
be used frequently. 

Lemma 3.4. [HI Chapter 5.8, Theorem 1] Let 1 < p < oo and C 1" 

a bounded Lipschitz domain. Then there is a constant C = C(p, n,Q) > 
such that for any v € W 1,P {^1) 

IMIp ^ c||Vi>|| p + ll^lli- 

We finally recall a "uniform version" of Gronwall's inequality. 

Lemma 3.5. |39|. Chapter III, Lemma 1.1] Let f,g,h be three nonnegative 
locally integrable functions on (io,oo) such that f'(t) is locally integrable on 
(to,oo), and which satisfy 

f'(t)<g(t)f(t)+h(t), t>to, 

rt+r ft+r ft+r 

I f(s)ds < ai, / g{s)ds < 02, / h(s)ds < 03, t > to, 

Jt Jt Jt 

where r, ai are positive constants. Then 

fit + r) < ( — + 03) exp(a 2 ), t > t . 



For the case p = 2 the initial data (vP,c°) has to satisfy (u°,c°) G 
W ' a (Q) x W 1,2 (Q) for the local strong solution (u, c) to be well-defined 
by virtue of Theorem 11.11 For simplicity we impose the following additional 
initial regularity condition on the concentrations for the rest of this article: 

Note that the L°°-constraint is not a restriction from the physical point of 
view, since concentrations of dissolved species are naturally finite. From 
the mathematical point of view, it is still possible to show that all results 
concerning global well-posedness and asymptotics in two dimensions pre- 
sented in this work essentially remain valid, if we allow for initial data 
c°) G Ll(n) x L 2 (n) with c° > 0. This will be explained in Remark[l3] 
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Definition 3.6 (Global solution). Let Q C R n , n > 2. The function 
(«, c) : [0,oo) — > L^(Q) x L 2 (Q) with q > is called global solution to 
(P), if it satisfies (P) almost everywhere and if for all T < oo it holds true 
that 

(u, c) G VF 1 ' 2 ^, T; L 2 (0) x L 2 (0)) n L 2 (0, T; V{A S ) x VF 2 ' 2 (f))). 

Let from now on Q, C M 2 be bounded and smooth. For the proof of 
Theorem 11.21 we aim for an estimate of type 

\\(u,c)\\ mTt2 <C(1 + T) 

with Et,2 the space of solutions defined right before (I2.24p . This requires 
appropriate estimates of the local solution (it, c). Applying energy meth- 
ods we will derive several uniform- in-time bounds for quantities like ||it||2> 
||c||2, etc., e.g. ||tt(i)||2 < C for t > and C > not depending on t. A 
justification for the use of energy methods, i.e. the formal multiplication 
of certain equations of (P) with quantities like u, c and integration over f2, 
can be given by the fact that (u, c) satisfies (P) almost everywhere and the 
integrals in the calculations being finite. Note at this stage that since we 
deal with strong solutions these boundedness relations are merely valid in 
an almost everywhere sense. However, as we aim at controlling the norm in 
Ey^, this is sufficient for our purpose. 

Let from now on (u, c) G Ej^ denote the unique local strong solution 
from Theorem 11.11 for initial data 

(u°,c°) g W$(n) x (w rl ' 2 (0)nL oo (0)). 

For (u,c) G L 2 (0) x L 2 (Q) with a > we define the energy functional 
1 N f 

E(u,c) :=-||n|| 2 + ^ / a log + |||V</>|| 2 + y ||0|| 2 ^, (3.1) 
i=l ^ n 

where is the solution to 

N 

— eAcj) = ^2 z j c ji d v (j) + t4> = (3.2) 

3=1 

Recall that £ is assumed to be a time-independent function on the surface 
dQ. Note also that the term Ci log Cidx is finite for Cj G L 2 (Q) with q > 0, 
since 

xlogx<(°' (3.3) 
In the following we set (xlogx)u = o := 0. 
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Let 

T ma x be the maximal time of existence of solutions (u, c). For t E 
[0,T max ) we define 

V(t) :=E(u(t),c(t)) 

N 



^ik*)ii!+E / c i {t)\o g c i {t)dx+*-\\vm\i+^ tt 



(3.4) 



The term Jo Cj log Cidx is bounded from below, since xlogx > —1, say, for 
x > 0. As the remaining terms in V are nonnegative, V is bounded from 
below. 

Lemma 3.7. Suppose (u,c) is the local strong solution to (P) from Theo- 
rem \l.l\ 

(a) We have 

^ / TT^I^WI 2 ^ e L} oc ([0,T max )), 
Jn DiCi{t) 

where ji = —D^SJci — DiZiCiVcf) denotes the flux of species i. 

(b) The function V from |3.^| ) is a Lyapunov functional for system (P). 
More precisely, it holds true that 

H N f 1 

_ m= -iiv„( () iii-y:y iiI ^u i ( ( )i^<o (3.5) 



for almost all t £ (0, T„ 



In the following proof we suppress all arguments of functions under inte- 
grals for better readability, whenever this does not lead to misunderstand- 
ings. 

Proof. ( of Lemma 3.1). In view of the derivation of (|3.5p and the mere non- 
negativity of Cj we have to first give a meaning to the term J n jj-^-\ji\ 2 dx, 
since it might not be well-defined. As will become clear in the proof below, 
this integral is connected to the term ^ J n a log Cidx, if it exists. We cannot 
just interchange differentiation and integration and then differentiate under 
the integral, since the derivative of the function x \— > x log x tends to oo for 
x — > 0. So it is not obvious that the term q log Cidx is differentiable with 
respect to time at all, since we only know Cj > 0. This is the reason why we 
first consider an approximation of this term. For 1 < i < N and 5 > let 

9i,s(t) ■= / (c; + 5) log(cj + 5)dx, t e [0, T max ). 
Jn 

Clearly is differentiable with respect to time and we have 

-Jl9i,s(t) = [ d t Ci log(cj + 5)dx + / dtadx. (3.6) 
d* Jn Jn 
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The last term is zero by (12.610 : with (jl.29H . (jl.31|) . integration by parts, the 
notation in (11.350 and straight forward manipulations we obtain 

;9i ) s( t ) = ~ [ ( div J i) lo g( c « + d) dx = [ Ji • V log (Q + 5)dx 
Jn Jn 

:Vq • Jidx 



dr 



Ci + 5 



[ - , 1 — T-(DiVci + DiZiCiV(t>) ■ (- AVcj - DiZiCiV <p)dx 
Jn DM + o) 

+ Zi °] ; V0 • (AVq + DiZidVcf) - Ciu)dx 
Jn c i + <5 

+ / — — — -u ■ (ziCiVcp + VcAdx 
Jn c i + o 



ji\ 2 dx — zi \ — l —V(f)-Jidx 



Di{ci + 8) JnCi + S 

f Ci f Ci 

+ Zi / — —^CiU ■ Vipdx + / — — rVcj • udx, t £ (0, T max ). 
JnCi + d J n Ci + b 

(3.7) 

Integrating equation (|3.7[) form to t £ (0,T max ) yields, after rearrange- 
ment, 

\ji\ 2 dx = gi,g(0) ~ 9i,s(t) - Zi [ [ °] • Jidxds 



>/!! 



Di{Ci + d) 7 Ci + 6 



+ z i / — C ——^iU ■ Vc/)dxds + / / — rVcj • udxds. 
Jo JnCi + o J JnCi + d 

(3.8) 

From this equation we are going to conclude assertion (a). The monotone 
convergence theorem implies that 

\ji\ 2 dx = lim / / 1 \ji\ 2 dx, 



'o Jn D iCi <5^o J J n Di(a + 5) 

however this limit might be infinite. To exclude this case we consider lim^oo 
of the right-hand side of (|3,8p . Relation (|3.3p and the fact that x\ogx > — 1, 
say, implies 

|(q + J) log( Ci + < l + ( Ci + 5) 2 , 

so there is an integrable majorant due to Cj(i) £ L 2 (f}), thus by continuity 
of x i — ^ x log x, x > 0, 

rimg i(5 (£)= / CilogCidx, te[0,T max ). 

For the remaining terms, note that Cj/(cj + 5) < 1, so by 

/ / \V(j) ■ Vci\dxds < ||V0]|i3(o,t;i2)||Vci||£,a C o 5t;i 2) < oo, 
io in 



o Jn 

t 



o Jn 
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Ci\V4>\ 2 dxds < ||ci|| L 2( O)t ._ L 2 ) ||V0|| i 4(o )t ; i 4) < oo, 
CiU- V(f>\dxds < 1 1 Ci\\ L 2(o,t,L2) I (0,t;i 4 ) 1 1 11^4(0,^4) < oo, 
f [ |VCj • U\dxds < 1 1 Vcj 1 1 x,2 ( ,*;X,3) 1 1 « 1 1 X,3 (0,*;X,2 ) < oo, 

io Jn 

the premises of the dominated convergence theorem are satisfied. Because 
of Ci/(ci + S) — > 1 as S — > taking the limit on both sides of equation f|3.8j) 
gives 

/ / 77— \ji\ 2 dx = I c° log c°cfe- I a(t) log Ci(t)dx 
Jo Jn UiCi Jn Jn 

/ V</> • Jidxds + Zi / Qu • V(f>dxds < 00, 
Vn io Jn 

(3.9) 

where we used the fact that due to u G the term Vq • ttcfcc = 0. 

Claim (a) is proven. 

Relation (|3.9p implies that 1 1— > J n C{ log CjC&r is absolutely continuous and 
its a.e. derivative is given by 

A / closed, = - / / V4>-M* + z, f CU-V^x (3.10) 

di 7n Did Jn Jn 

for almost all t G (0,T max ); note that the right-hand side of (|3.10p is in 

£L(°> T ma*) duet ° 

In two dimensions the energy equality for weak solutions of the Navier- 
Stokes equations is valid, cf. [371 Theorem V.l.4.2]. Since strong solutions 
are also weak solutions, we have 



1 A N 

II 2 

2d7" M " 2 



N 

V^lll - y^ z » / CiV(j)-udx, t E (0, T max ). (3.11) 



Recalling the definition of V in (|3.4p with the help of summation of (|3.10p 
over i, (|3.1ip and interchanging differentiation and integration, we obtain 



d N f 1 N f 

—v(t) = -\\vug - v / -— ijipds-y;^ / v<p-Jidx 

+ e / V0 • V(j)tdx + et 4>4> t dx 
Jn Jan 



(3.12) 



for almost all t G (0, T max ), where we write 0j := 5*0. Note that the 
respective last terms of (|3.10p and (|3.1ip cancel. Observe also that due to 
the time independence of £ differentiating the boundary condition of (|3.2p 
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in time yields d u (frt + Tcpt = 0. Hence using integration by parts, taking into 
account (ll.29|) , (jl.31j) and the time- independence of £ gives 

2_] Zi / V(f> • Jidx = / 4>dt 2_] z i°i ) dx = —e / <f)A(f) t dx 
i= i J » Jo. \i=i / Jn 

(3.13) 



e I V<fi • ^4>tdx — e / 4>d v (\>tdx 
Jn Jon 

e I V(p ■ Vc/)tdx + et I cf)(j)tdx. 
Jn Jan 



Thus, collecting f)3. 12|) and (|3.13[) . we see 

H N f 1 

±V(t) = -\\Vu\\l-Y, / j^\j t \ 2 dx<0 
dt f-j* J n Did 

for almost all t £ (0,T max ). □ 



Remark 9. Observe that due to V € W^ G ' C (0, T max ) the function 1/ is ab- 
solutely continuous, so for almost all t £ (0,T max ) it holds true that 

v(t) = v(o) - £ f ||v«wnl + J n D^)^ 2dx ) ds - ( 3 - 14 ) 

Hence V is non-increasing. 

Lemma 3.8. There is a constant C > 0, depending only on the initial data, 
such that for almost all t £ (0, T max ) 

|K*)|| 2 + ||ci(t)]ogCi(*)||i < c. 
Proof. From ()3.14j) we have V(t) < V(0), hence 

K*)ll2+ / Ci(t)\oga(t)dx < C , 
Jn 

where Cq > only depends on the initial data. From CjlogCj > — 1 we 
deduce 

/ \ci(t) log Ci(t)\dx < / Q logCidx + 2\Q\ < C' , 
Jn Jn 
which shows the claim. □ 

With Lemma [3721 the boundedness of the term ||cjlogCj||i enables us to 
show that the L 2 (fi)-norm of concentrations Cj remains uniformly bounded 
as long as the solution exists. 

Lemma 3.9. There is a constant C > 0, depending only on the initial data, 
such that for almost all t £ (0, T max ) 

<C and f \\Vci(s)\\lds < C(l + t), i = l,...,N. 
Jo 
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Proof. We multiply equation (|1.29p by Cj and integrate over ft by parts, 
making use of (ll,3ip . to get 



Id 
2d? 



— — ^ 1 1 Ci 1 1 1 + Dj||Vcj||2 = — DiZi I qV <p ■ V Cidx + \ CiU-Vcidx. 



Since 



/ CiU-Vcidx = - [ u-V(c 2 )dx = 

due tone Ll(Q) and V(cf) G L 2 (Q), we are left with 

\ci\\l + Di\\Vci\\l < C / \ciV(f)- Vci\dx. (3.15) 
Jn 



1 d_ 

2 dt 1 



In order to use a Gronwall argument we have to estimate Ii := j n |cjV<^> • 
Vci\dx appropriately. For convenience we reproduce the computation from 
|10j correcting also a slight mistake. We use the same splitting into <p\ and 
4>2 for the potential <f> as in (|2.ip and (|2.2p . Let 1 < q, r < oo be such that 
1/2 = l/<? + 1/r. Then 

< ||Vci|| 2 ||ci|| r ||V^||g < ||Vci|| 2 ||ci|| r (||V^i|| 9 + \\V<h\\q)- (3-16) 

With W 2 ' 2 {Q) ^ W^(n) for any q G [l,oo), since ft C K 2 , due to 
Sobolev's embedding theorem and the help of Remark [2] we compute 

||V0i|| g < \\<f>i\\w^ < C\\<h\\ww < C'U\\ w y 2{my (3.17) 

hence ||V0i|| g < C independent of time. Now let 5 G (1,2) and q := 
then by Sobolev's embedding theorem we have W 2 ' S (Q) W 1,q {Q,). So 
again taking advantage of Remark [2] we have 

N 

\\V<h\\q < \\H\W^ < C\\<h\\w*.> < C 'J2 IM*" ( 3 - 1§ ) 

J'=l 

Combining (f3TT6l) - ([3TT8]) yields 

h < C\\Vci\\ 2 \\ci\\ r f 1 + J] H c jH* j ' ( 3 - 19 ) 

where the constant C does not depend on t or q. 

Applying Lemma 13.11 and Lemma 13.21 and the fact that ||cj||i = mj and 
\\ci log Q ||i < C due to Corollary Q] and Lemma [3 .8^ we see that for any p > 1 
and e > there is a constant C > such that 

||Ci||p 5: e ||^ c «|l2 ^ + C) 
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where C depends on mj and the parameters indicated in Lemmas 13.11 and 
Thus we may estimate I{ by 



h < C||Vc i || 2 (£||Vci||2~ 1/r ' + C) [e\\ E l|V Ci ||2~ 1/5 + C 

< C (l + g || VcJ^) + £ g ||Vc,||^ + 1 

(3.20) 

for e < 1. With the help of 

Af N 
i=l j=l 



N 

,3— 1/5 
-ill 2 

i=l 



< cEIIVq 



and Young's inequality we obtain 

N , N 



X>< c( e E(HVQ|l2" 1/r " 1/5 + l|Vc,||I) + l) . (3.21) 

i=l ^ i=l 



From n = 2 and the definition of r and 5 we see 1/r + 1/5 = 1. So, writing 
d' := min{Dj} and choosing e > small enough, we infer from (13.211) 



d> N 



oEH Vc ^l2 + C ( 3 - 22 ) 



2 

i=i i=i 
Summing (|3. 15)) over i and using estimate (I3.22H yields 

N N 



1 ^ - d, Eii Vc *ii2 + c - ( 3 - 23 ) 

i=l i=l 



Poincare's inequality (cf. Lemma [3. 4|) implies that there is a constant C > 
such that 

-II Vci|| 2 < -C||ci|| 2 + C||ci||i < -C||q|| 2 + C. 
Plugging this into ()3.23|) results in 

, N TV 

i=l i=l 

With Gronwall's inequality we deduce that there is a C > such that 

\\ci(t)\\ 2 < C, t€[0,T max ). 
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Integrating (|3.23j) from to t gives 



r||Vcj(s)||lcfo<C"(l + i) 
Jo 



□ 



Remark 10. If we integrate (|3.23p only over (t, t + r) for t < T max and 
r > small enough such that the integration makes sense, we obtain 

\\Vc t (s)\\ 2 2 ds<C(r) (3.24) 

independently of t. This will prove important for the application of the 
uniform Gronwall inequality in Lemma 13.51 



Remark 11. In three dimensions it is also possible to estimate the term i, 
along the lines of the proof of Lemma [3.91 using Corollary [T] and Lemmas 13. 1| 
13.21 and 13.81 However performing the same calculations one observes that in 
this case 

i + ±<i, 

r o 

so ij cannot be absorbed into ||Vcj|||, cf. [10j . 

Lemma 3.10. There is a constant C > 0, depending only on the initial 
data, such that for almost all t £ (0,T max ) 

+ \\<Kt)\\w,» < c 



Proof. Note that once the uniform boundedness of ||cj||4 for t > is obtained, 
by Sobolev's embedding theorem and Remark [2] it follows that 



W i >00 ^ ^\m\ w 



< C" E NU + 1 \< C " & € (0,T max )). 



Let k e N and t > 0. We multiply (fL29]) by c 2k "\ integrate over and by 
parts and use again the no-flux condition (|1.3ip to obtain 

= (2k -1) J (-D iC l k - 2 \V Ci \ 2 - DiZicf _1 VciV^ + cf ~ l uVc^ dx 
= -A^r=^ / \V(c k )\ 2 dx-D iZi ^^ ! c k V(c k )V<bdx 

(3.25) 

To make this formal calculation rigorous, it is important to note 
that (u, c) satisfies problem (P) almost everywhere and that Cj G 
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PF 2 1/2 (0,T;^ 2 (fi) ^ L*(0,T;L*(n)) for any q < oo, thus ||t£|| 2 = \\ci\\% k 
is well-defined for almost all t £ (0,T max ). 

Due to the divergence free condition on u the term f Q u • V(cf k )dx is zero. 
If we denote 

v k := ct (3.26) 

we can rewrite (13, 25ft as 

Id 2k — 1 2k — 1 f 

-^Wh-A- ^— H Vu fclli- A*— — J v k Vv k -Vcj>dx. (3.27) 

The Holder inequality yields 

1 H 2A- — 1 2fr — 1 

2jfc^ll«fcHl < -A^^IIV^Ill + C— ^- ||«*||3||V«*|| a ||V^|| 6 . (3.28) 



The uniform boundedness of 1 1 1 1 2 due to Lemma 13.91 and the Sobolev em- 
bedding W 2 > 2 (n) ^ W lfi {Q) imply 

\\v^ 6 <M w ^<c\^2\\ Cj \\ 2 + ?j <c. 

Applying Lemma 13.31 to || u^. || 3 shows that 

hkh < C(ll V-l; fc|l2 /3 ||^A:|li /3 + \\Vk\\l), 

where C > is independent of v k and t. Hence we estimate (I3.28|) further 
by 

^Nlli < l|v«*lli + c||v^|| 2 (||v^||f \\v k \\\ /3 + Kiio. 

To obtain the lemma we set k = 2; note that ||i>2(i)lli = Il c «(^)ll2 — so 
~ t \\Ml < " y l|V^||I + C'UV^IhdlV^llf + 1) 
<-^\\Vv 2 \\t + C", 

where we made use of Young's inequality in the last step. Applying 
Poincare's inequality (Lemma I3.4|) and taking into account the bounded- 
ness of 11^2 ||i yields constants d',C > 0, such that 

^IMIl < -d'\\v 2 \\l + c. 

Finally, from Gronwall's inequality \\v2W2 = IIqIII is uniformly bounded for 
t > 0, since ||cj(0)||4 = 1 1 c° 1 1 4 < 00. □ 



We are now in position to obtain suitable estimates on the velocity field u. 
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Lemma 3.11. There is a constant C > 0, depending only on the initial 
data, such that for almost all t € (0,T max ) 



/ ||V«(s)|||ds < C(l + 
Jo 



<0 

Proof. Recall the energy equality from (|3.11j) . i.e. 

N 



1 d o 2 

-~r \\u U + \\vU U 
2di" 112 11 112 



/ V ZicjVift • udx. (3.29) 
Jn t=1 



The right-hand side can be estimated with the uniform boundedness of \\u\\2, 
\\c\\2 and ||V</>||oo, cf. Lemmas 13.81 13.91 and 13.101 as follows: 

dV(j)-udx < \\ci\\ 2 \\V(f)\\oo\\u\\2 < C (t £ (0,T max )). 

So integrating (|3.29[) from to t in time gives 

+ f l|V«(a)||Id S < i||n°|| 2 + Ct < C'(l + t), 
which proves the claim. □ 



Remark 12. Note that if we choose t and r in the same way as in Remark [TU1 
and integrate (|3.29|) over (t,t + r), this results in 

J + ||Vit(*)||!ds < C'(r), (3.30) 
where C only depends on r, but not on t. 

We also need estimates on higher derivatives of u. 

Lemma 3.12. There is a constant C > 0, depending only on the initial 
data, such that for almost all t £ (0,T max ) 

||V«(t)||i < C and [ \\A s u(s)\\ 2 2 ds <C{l + t) (3.31) 

Jo 

where As denotes the Stokes operator defined in Remark{J\ 



Proof. Multiplying (|1.28j) by A$u and integrating over Q yields 

1 d 

2dt 



Id 



Vu\\l + \\Asu\\% = - / (u ■ V)u ■ A s udx - / ZjCjVcf) ■ (A s u)dx 
Jn Jn i=1 

(3.32) 

where we used the well-known fact that 

u t . A S U4 X = - f « t AWx = f Vu, : Vu - [ , A - = ^l| V* 
Jn Jn J dn 2 at 

(3.33) 
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Working in two dimensions, the Gagliardo-Nirenberg inequality implies 

\\v\\a < C|M|2 /2 ||V?;||2 /2 , for v G W 1,2 (£l). (3.34) 

The Stokes operator in L^(f2), C W 1 a bounded smooth domain, is closed 
with G p(A s ), cf. [33 Theorem III.2.1.1]. Hence, we infer ||V 2 tt|| 2 < 
||n||^y2,2 < C||j45u||2 for u G T>{Ag). This in combination with (|3.34p . the 
uniform boundedness of 1 1 tt 1 1 2 by Lemma 13.81 an d Young's inequality yield 
the following estimate on the non-linearity. 



(u ■ V)tt • Asudx 



< \\u\\ A \\Vu\\i\\Asu\\ 2 
2 



< c|M|i /2 

< ±\\A s u\\ 2 2 + C>\\Vu\\i 



i\\ 2 \\V z u\\ l 2 /2 \\A s u\\ 2 



(3.35) 



Using the a priori estimates on Cj and V0 from Lemma 13.91 and 13.101 and 
Young's inequality the nonlinear electro-kinetic part can be treated by 

N 



Jn i=1 



< C7||V0||oo||c|| 2 ||^t*I|s 



<h\A s u\\l + C. 



(3.36) 



So estimating ()3.32j) with f)3.35|) and (|3.36p results in 
d 



At 



\Vuf 2 + \\A s uf 2 < C\\Vu\\\ + C. 



(3.37) 



Setting g{t) := C|| Vii||| this implies 



d 
di 



\Vu\\ 2 2 <g(t)\\Vu\\l + C (fe (0,T max )). 



If we let f(t) = llVulU, and h{t) := C, with Remark 1121 all premises for the 
uniform Gronwall inequality Lemma 13.51 are fullfilled, and hence 

||Vu||l<C" (t G (0, T max )). 

Integrating ()3.37p over (0, t) therefore yields 



\\A s u{s)\\lds < C\\ + t). 



□ 



The estimation of higher derivatives of Cj turns out to pose certain difficul- 
ties since a calculation analogous to (|3.33p . creates, in general, non- vanishing 
and unpleasant boundary integrals. Indeed it gives 

/ d t Ci-Acidx = - / V(d t Ci) -Vcidx + / (d t Ci)(d v Ci)dx 
Jn Jo, Jan 

(d t Ci)(dyCi)dx. 



1 d „ 
Vc 

2 At 



2 

i 112 
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By this fact, we prefer to introduce new variables instead of working with 
Ci directly. In |10j a problem similar to (jl.29j) - (jl.31jl is transformed to a 
problem for the new variable 

:= Ciexp(zi<j)). (3.38) 

We will show that £ possesses exactly the regularity as desired for c in 
Theorem ll.il i.e. 

C G W x ' 2 (O)(0, T; L 2 (fl)) n L 2 (0, T; W 2 '\fl)), 

cf. Lemma [3.151 For, we plug Cj = Q exp(—Zi(p) into (|1.29p - (|1.3ip and obtain 
the following nonlinear heat equation subject to homogeneous Neumann 
boundary conditions for Q: 

dtd ~ DiAQ = -DiZiV<j) ■ VCi + Zid4>t + ZiQu ■ V0 - u ■ in Q, 

(3.39) 

d u & = on dn, (3.40) 

for i = 1, . . . , N. Those homogeneous Neumann boundary conditions enable 
us to apply the corresponding relation to ()3.33p for £, namely 

/ dtd ■ A(idx = - [ V(dtCi)-V(idx+ [ (dtCi)(d v (i)dx = -~\\V(i\\l 
Jn Jn Jn * a* 

(3.41) 

Once we have "good" estimates for higher derivatives of £ we will be able 
to infer the same behaviour for higher derivatives of c, since c and £ are 
connected via eft, which already has nice regularity properties. 

However, in order to obtain suitable estimates for V£ and A£ via the 
energy method, we need to gather more information on the terms on the 
right-hand side of (|3.39p . This will be done in two steps. 

First we shall prove the existence of a uniform L°°-bound in space and 
time for concentrations q (Lemma 13. 13|) . With this information at hand we 
gain sufficient information on the new term 4> t (Lemma 13. 14|) . 

Lemma 3.13. There is a constant C > 0, depending only on the initial 
data, such that for almost all t G (0,T max ) 

||Ci(f)||oo < C. 

The computation for the proof of Lemma 13.131 is essentially contained in 
|10j . for convenience we provide a proof here correcting a minor mistake. 
The idea is as follows: If we can show that there is a sequence pk — > oo 
with the property that esssup tg r 0iOO ) ||cj(£)|| Pfc < C independently of k, then 
not only is c;(t) G L°°(f2) for almost all t G [0, oo), but it also holds that 
Ci G L°°([0,oo) x Q). This method is sometimes referred to as "Moser 
iteration", cf. [4"3] . 
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Proof, (of Lemma \3.13\) . Recall the notation v k = for k E N. Taking into 
account the L°°-bound for V0, we improve the estimation of (I3.27|) given in 



(pT28l) and get 



^^Iklll < -^IIV^HI + C||t;*|| 2 ||V«*|| a . 



Thanks to Lemma 13.31 and Young's inequality, for every e > there is a 
C{e) > such that 

^ t \W\\l - f liv^lli + C\\Vv k U\\Vv k \\T \\v k \\\ /2 + |K||i) 

^-^llV^lll + ellV^lll + C^H^II?. 

If we choose e = & , we get 

^HWi < -§H V ^lli + C'(k)\\v k \\l (3.42) 

where C = C'{k) can be controlled by a polynomial in k. Poincare's in- 
equality (Lemma 13. 4D implies 

-llVffeHa < -7HKII2 + Ikfelli 

with a constant Cp > max{l, without loss of generality. Plugging this 
into (|3.42p and multiplying with 2k gives 

^IKIll < -^-llwfclll + (2kC'(k) + Di)\\v k \\l (3.43) 
Considering a fixed t £ (0,T max ), we define 

r / /■ , 

5fc := max < ||cj(0)||oo,esssup I / c i (s,x)dx 

0<s<t \Jq 

Note that we already proved that S±, S 2 , S4 can be bounded independently 
oft. 

With the abbreviations C\ = Di/Cp and C 2 {k) = 2kC'(k)+Di we deduce 
from equation (|3.43p that 

^ t \\v k (t)g < -Ci||t; fc (t)||! + C 2 {k)Sf, 0<t<T. 
Multiplying with e Clt yields 

dt 

By a standard Gronwall argument we therefore can achieve that 



-\\v k {t)\\l<C 2 {k)Sl k e c ^. 



k(*)lli < K(o)||l + ^sf < m\\4\\ 2 * + 

< I in, + 9m S f < 2^)5*, < t < r max , 



(3.44) 
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since we assumed that Cp > Therefore 

\ 1 / 2k i/o, /ir„(h\\ l / 2k 



a 



c 

n 



i(t,x) 2k dx^j = (\\v k (t)\\lf 2k < (^^) S k . (3.45) 



We can find constants C, m > such that 



2C ^ k ) - 2Cp ( l + 2k 9M\ 



Collecting (I3.44| )-( |3.46p . we obtain the following recursion relation for Sk- 

5 2fc <ma X ||| Cj (0)|U^^-j S fc j = (-^-j 5fc (3.47) 
< C 1/2k k m/2k S k . 

Thus from the fact that £4 < 00, we deduce that 5g < 00 and so on. Note 
that this recursion is valid for all < t < T max , and the constant C does 
neither depend on t nor on k. 

If we set k = 2^, fi = 0, 1, 2 . . ., relation (|3.47j) implies 

S 2 , + i<C 1 / 2 " +1 2 m ^ 2li+1 S 2 ,, ix = 0,1,2,.... 

Since 5i < 00 due to Corollary [T] and the series ^2~^ 1+1 ' and ^/x/2 M+1 
converge, the sequence (SV)^ is bounded. Note that this means that 

ess sup || c(i) || 20 < 5V < C 
independently of fi. We therefore proved ess supq<^ < T max 1 1 ^Wl loo <c. □ 
Let us consider <f> t = dt4>; for almost all t E (0,T max ) it holds true that 

N N 

—eA(j} t = z i&tCi = — z i °-i v x e (3.48) 

i=l i=l 

d v <k + T (j) t = 0, x E dn (3.49) 

by interchanging <9 t and A and taking advantage of f|1.29j) and the time- 
independence of £. 

Lemma 3.14. There is a constant C > 0, depending only on the initial 
data, such that for almost all t E (0,T max ) 

\\Mt)h<c. 

Proof. The proof follows by a duality argument. Let 

Bh := -eAh, h E V(B) := {v E W 2 ' 2 {n) : d v v + tv = 0}. 

Remark [2] then implies that B : D{B) — > L 2 (Q) is isomorphic. By duality 
also B* : L 2 (f2) — > V(B)* is isomorphic. By the self-adjointness of B and 
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since V{B) is dense in L 2 (Vt) we further can regard above B* as the unique 
extension of B, hence we write 

B : L 2 (Q) -> V(B)*. 

According to Theorem ll.il we have 

div Ji G L 2 (fi) «^ V(B)*. 

We now intend to show that the norm of each div Jj in D(B)* can be esti- 
mated uniformly in t G (0,T max ). Thanks to (|1.31[) we can calculate 



(div Ji,h) V ( B y = / hd\\Jidx= I Vh-Jidx 
Jn Jn 

= / V/i • (— DiVci — DiZiCiVcj) + Ciu)dx 
Jn 



(3.50) 



for all h G T>(B) and for almost all t G (0,T max ). 

We need to estimate the terms on the right-hand side of (|3.50p . Applying 
integration by parts and Holder's inequality we compute 



/ 

Jn 



Vcidx 


< 


I C{Ahdx 


+ T 


/ Cihdx 






Jn 




Jon 



< 



Hlr>(B)IMl2 + T||ci|| 2) 9n||^||2,an 



< C(IM|2 + ||Cj||oo) 

<C\\h\\ v{B) (t G (0, T max )), 

where we made use of the uniform boundedness of ||cj||2 and ||cj||oo) cf. 
Lemmas 13.91 and [3.131 With the uniform bound on || V</>||oo, cf. Lemma r3.101 
we deduce 



CjV0 • Vhdx 



< Ik,; 



WU\\h\\wi,i<c\\h\\ v(B) 



for the second term. Due to Lemma 13.81 the quantity \\u\\2 is uniformly 
bounded in time, so the last term can be estimated by 



au ■ Vhdx 



< lie,- 



\u\\ 2 \\Vh\\ 2 <C\\h\\ v(B) . 



In total we obtain 

N 

Uth < £ div J ih <CJ2 II div J iWv{By <C (te (0,T max )). 



N 



i=l 



i=l 



□ 



Before treating problem (|3.39p - (|3.40p let us note the following lemma on 
the regularity of £. 

Lemma 3.15. Let c G E^ 2 be as in Theorem \1.1\ and ( be defined as in 
fPSj) . Then it holds true that ( G E^ 2 . 
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Proof. This is an immediate consequence of 4>, V</>, Cj and, hence, A0 being 
uniformly bounded in L°°(i7) and <^ being uniformly bounded in L 2 (f2) for 
any time. □ 

Lemma 3.16. There is a constant C > 0, depending only on the initial 
data, such that for almost all t £ (0,T max ) 

\\Vci(t)\\l<C and [ \\A Ci {s)\\lds <C(l + t). 

Jo 

Proof. We first derive the corresponding estimate for the transformed vari- 
able d, defined in We multiply (pTBUD with -AQ and get with (ET4"T]) 

1 d_ 

2ch 



— Zi QiU ■ VcfiAQdx + / u ■ VQA(idx 



(3.51) 



Note that the L°°-norm of Q is bounded independently of time, since both 
|| c|| oo and ||<^||oo are bounded independently of time, cf. Lemmas 13.101 and 
13.131 Estimating the integrals on the right-hand side, apart from the last 
one, with Holder's inequality, Lemmas l3.8U3.12l 13. 14| and Young's inequality 
gives 

V0 • VCiAQdx < ||V0||oo||ACi||2||VCi||2 < exIIACiHl + C(ei)||VCi||!, 

(ifaMidx < HCillooll^tlbllACilb < e2||ACi||| + C(e 2 ), 

Qu ■ VcfrA&dx < ||Ci||oo||V0||oo||«||2||ACi|| 2 < esllACilH + C(e 3 ). 

(3.52) 

For the remaining integral f„ u-VQAQidx recall that if v £ W 2,2 (Q) satisfies 
homogeneous Neumann boundary conditions, one has ||V 2 i> H2 < C||Av||2. 
So Holder's inequality and the Gagliardo-Nirenberg inequality (|3.34p give 



u-VQAddx < ||it||4||VCi||4||ACi||2 

1/2 ||Vn||f HVOIlf "^ l|1/2 
<£ 4 ||ACi||i + C( £4 )||VCi|||, 



< Cllulli^HVulli^llVCilla^llV^iliniACilh < C"||VCi||2 /2 ||ACi||2 /2 

(3.53) 



where we also applied Lemmas 13.81 and 13.121 

Collecting f|3.52j) and f|3.53j) and choosing £j small enough we obtain from 



dt 



||VCi||l + A||ACi||l < C\\VQ\\ 2 2 + C. 



(3.54) 
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It holds that VQ = Vq exp(zi<f>) + ZiCi exp(zj0)V^>. Since 4> and are 
uniformly bounded in time and space, V£ has the same growth as q in 
(pOIjh namely J f <+r ||VCi(s)|||ds < C(r), t > 0. Letting f(t) = ||VCi|||, 
#(£) = /i(t) = C in the formulation of Lemma 13.51 we find 

||VCi|||<C" (t G (0, T max )). (3.55) 

Integrating (|3.54p in time finally gives 

f \\ACi(s)\\ 2 2 ds<C'(l + t) (3.56) 
Jo 

for some C > 0. 
Since 

Vq = VCiexp(-2: i 0) - Zi(iV(pexp(-Zi(p) and 

Aq = Adexp(-Zi(j)) - 2ziVd ■ V^>exp(-^0) 

- ZjCiA(/)exp(-Zj^) + Zj 2 Ci|V(/>| 2 exp(-2; i 0) 

the estimates (|3.55|) and ([3.56P are also valid for Vq and Aq by the uniform 
L°°-boundedness of (j>, V</> and Q. □ 

The estimates on (u, c) and on <f> derived so far enable us now to conclude 
the statement of Theorem 11.21 

Proof, (of Theorem\TE). From Lemmas EES EES I3TT21 and [37T61 it follows 
that 

II U IIl 2 (0,T;D(A s )) + W C \\l 2 {0,T;W 2 - 2 {Q)) - + T )- (3.57) 

We use this estimate to show the boundedness of the nonlinearities of (P) 
in order to prove the boundedness of <9j(u, c) in L 2 (0,T;L 2 ). We have: 

\\ c i c j\\h(Li) + ll Vc * • V0|||2 (L 2) < CT. 

For the treatment of the term u ■ Vq we will, as in the proof of Lemma [3.161 
work with the transformed variable Q. The reason lies in the homogeneous 
Neumann conditions for Q which allows for the estimate ||V 2 (j||2 < C|| A£j||2 
when using the Gagliardo-Nirenberg inequality. Since 

u • Vq = u • VCi exp(— Zi4>) — ZiQu • V</>exp(— Zi4>), 

the estimate 

f T 

\\ u ■ v Ci|lL2 (L 2) < / ||u|| 2 ||V0||idt 



< C [ ||u|| 2 ||Vu||2||VCi||2||ACi||2* < C(l + T), 

J 



is sufficient to show \\u ■ Vq|| 2 2 ( i2 ) < C{1 + T), since £j, and V0 are 
uniformly bounded in L°° , cf . Lemmas 13.101 and 13.131 and u is uniformly 
bounded in L 2 , cf. Lemma 13.81 We infer 

||^ Cl || i2(i2) <C(1 + T). (3.58) 
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We proceed in the same way for the Navier-Stokes part, i.e. 

r T 

\\P(u ■ V)ti||^ 2{L2) <C \\u\\l\\Vu\\ldt 

r T 

<C ||n||2||Vn||^||V 2 n||2(it 
Jo 

< C" ( \\u\\2\\Vu\\ 2 2 \\A s u\\ 2 dt < C""(l + T), 
J o 

N 

|| Y, ZiCiVHh^) <CT, 

i=l 

where we used the fact that due to E p(A$) we have ||V 2 u||2 < C||As«||2 
for u E V(Ag), hence also 



| L2(L2) <C(1 + T). (3.59) 

In conclusion it follows from (|337|) - (|339| ) that \\(u, c)||e Tj2 < C(l + T) for 
finite T. So (u, c) can be extended to a global strong solution. It is unique 
by the local uniqueness from Theorem ll.il □ 

Remark 13. It is possible to show that a global existence result in the spirit 
of Theorem 11.21 remains still valid if we only impose (u , c°) E L 2 (Q) x L 2 (Q) 
with c° > for the initial data. This will be published in the doctoral thesis 
of the second author. We only briefly indicate a sketch of the proof. 

1. Existence and uniqueness of local-in-time weak solutions (u,c): 

The idea is to define a map L which sends c E C([0, T]; L 2 (Q)) n 
L 2 (0,T; W 1 ' 2 ^)) with c(0) = c° to the unique (weak) solution c E 
C(0, T; L 2 (n)) n L 2 (0, T; W 1 ' 2 ^)) of 

QjCj + div(-L>jVcj + Ciu) = DiZi div(cjV^), i = 1,. . . ,NA 

d u Ci\dQ = -(zic~id u (j>)\dn, i = l,...,N, > (3.60) 

c;| t=0 = c°, i = l,...,N, J 

i = 1, . . . ,N, where 4> is the unique solution to 

N 

— eAcj) = z jCj> in n, d v (j) + T(j) = on 9fi (3.61) 
i=i 

and it is the unique (weak) solution to 

d t u + A s u + P(u-V)u = -P^Y z i^ V ^i u(0) = u°. (3.62) 

The unique weak solutions u and c exist from well-known results on weak 
solutions for Navier-Stokes equations, see e.g. [37], and for parabolic equa- 
tions, see e.g. [27], whereas 4> is well-defined by Remark [2) By choosing the 
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time T small enough it is possible to show that L is a contraction in a closed 
subset of 

Z T = {ve C([0,T];L 2 (n)) n L 2 (0,T; W 1,2 (fi)),t;(0) = c }. 

Hence L has a unique fixed point, which represents a local-in-time weak 
solution. This idea has also been used in |35j . 

2. Usage of regularizing effects of Navier Stokes and parabolic equations: 

The goal is to infer that for some positive time the trajectory (u(-),c(-)) 
enters the regularity class for initial data considered in this work, namely 
Wo,ff ( n ) x (W 1,2 (fi) n L°°(fi)). This can be seen as follows. 

Since the weak solution (it, c) is contained in L 1 (0,T; Wq' 2 (Q) x W 1,2 (Q)), 
by Lebegue's differentiation theorem we obtain that almost every t G (0, T) 
is a Lebesgue point for (u, c). Now let to £ (0, T) be such a Lebesgue 
point. Then (u(t ),c(t )) G ^(O) x W 1 ' 2 ^) is well-defined. Thanks to 
Theorem 11.11 there is a unique local strong solution 

(u s ,c s ) G W^&T^LKn) xL 2 (O))nI 2 (0,T„;D(i s ) x ^ 2 > 2 (ft)) 

with (tt s (0), c s (0)) = (u(to), c(to)) and for some To < T. Note that strong so- 
lutions are also weak solutions, so by uniqueness of weak solutions they must 
coincide, i.e. (u(t +t),c(t + t)) = (u s (t) , c s (t)) for t G [0,T ). Therefore 
(u(t),c(t)) G WfflSl) x W 1 ' 2 ^) for all 4 G (to,t +T ). 

It remains to gain L°°-regularity for c. To this end, however, we can argue 
in the same fashion. From 

cGL^ ,*o + ?b;W 2 > 2 (n)), 

we see that c(h) G W 2 ' 2 (tt) L°°(fi) for a Lebesgue point ii G (t ,?b). 
Consequently, we have reached the initial situation of this paper. 



4. Steady states 

For the definition of steady states we only require that they solve (P s ) in 
the weak sense. 

Definition 4.1. (Steady state) Let Q C W 1 , n > 2. The function (u, c) G 
Wo' 2 (n) x W 1 ' 2 ^) is called a steady state of (P), if (u, c) is a weak solution 
to (P s ), i.e. if for all ((p,if>) G Cg^(fi) x C°°(Tl) the following is valid: 

N 

— / • V(pdx — / u(u • Vip)dx = \^ z j / CjVc/) • <pdx, (4.1) 
./n Jn - =l Jn 

/ (- AVcj - DiZiCiV(j) + Ciu) ■ Vipidx = 0, i = 1, . . . , N, (4.2) 
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where <fi G W 3 ' 2 (£l) satisfies 

N 

— eAcp = z j c j i n ^> 9 U 4> + T<fr = £ on dQ,. (4-3) 

7=1 

It turns out that steady states have in fact better regularity properties as 
required from Definition 14.11 

Lemma 4.2. Let (u°°,t^) G w££(Sl) x W 1 ' 2 ^) be a steady state of (P). 
Then (u^c 00 ) G W 2 > 2 (Q) x W 2 ' 2 \ti). 

We omit a proof here. It is an easy consequence of a bootstrap argument 
taking into account the semilinear structure and regularity properties of the 
Stokes operator and the Neumann Laplacian, cf., e.g., [16], [37], and [3D]. 

Let SI C I 2 be bounded and smooth. Before considering existence and 
uniqueness of steady states we give a characterization of weak solutions to 
problem (P s ). 

Lemma 4.3. The function (n 00 ^ 00 ) G Wq£(SI) x W l > 2 (£l), where cf > 
not identically zero, i = 1, . . . , N, is a weak solution to (P s ) if and only if 
the following conditions are fulfilled: 

(i) u°° = 0, 

(ii) C°° := cf exp(z i 0°°) = const, i = l,...,N, 

where is the corresponding electrical potential from |^.3[ ). In this sit- 
uation there is a constant r/o > such that cf > r)o, i = 1,...,N. The 
associated pressure tt°° satisfies 7r°° = ^jcf up to some constant. 

The proof of this lemma employs the functional V from (I3.4p . In terms 
of the transformed variable £ from (|3.38p . the dissipation rate 4;V(t) can 
be rewritten as 

-V(t) = -\\Vu(t)\\ 2 - T / ^e X p(-2^0(t)|VC,| 2 dx, (4.4) 

whenever (u, c) is a global strong solution from Theorem 11.21 

Proof, (of Lemma \4-3^ . Suppose conditions (i) and (ii) are satisfied. Since 
u°° = and 

= VC = (Vc°° + ZicfV^) exp(^°°), (4.5) 

we infer that (14, 2p is satisfied. From (14. 5h it is now evident that the term 
— Yl,j ZjcfVc/) 00 = V Ylj c j° is a gradient field, so also (]4.1j) is valid. From 
this the last assertion of this lemma concerning the pressure follows, too. 
Observe that (f = cf exp(zi</>°°) > 0, because q > and Cj ^ 0. The fact 
that G IF 3,2 (0) BUC{p) due to Sobolov's embedding theorem shows 
that there is a constant % > such that cf > r]o, i = 1, . . . , N. 
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Now suppose (u 00 ^ 00 ) is a steady state to (P). From Lemma 14.21 the 
function (u s (t), c s (t)) := (u 00 ^ 00 ), t > 0, is a global strong solution to 
(P) in the sense of Definition 13.61 For this stationary solution we have 
V(t) = E(u°°,c°°), hence ^V(t) = 0. On the other hand, with flU} we 
deduce Vii°° = and VC°° = 0, which results in (f° = const and u°° = 
because of u°°\qq = 0. □ 

In order to make use of the global strong solutions (u, c) from Theorem ll.2l 
it is crucial to know that the dissipation terms in (|4.4p become small as time 
tends to infinity. 

Lemma 4.4. Let (u,c) be a global strong solution to (P) from Theorem \1.2\ 

and Q = c% exp(zj</>) be defined as above. Then 

lim (||V^)||I + ||VC^)||I)=0. 

Proof. From (|4.4p and the uniform boundedness of ||</>||oo an d ll c lloo> cf. Lem- 
mas 13.101 and 13. 13|, we have 

i N 1 

±V(t) < -\\VU\\1 - J2 7 9 TTH IIV6H1 

at ^ l|ci||oo|| exp(-2^0)||oo 



<-c(\\Vuf 2 + J2\\Vb\\£j 



The functional V is bounded from below, so using the fact that it is abso- 
lutely continuous, cf. Remark [9j we compute 



L < V(t) < V(0) ~ c J q ^ll Vn ll2 + YlW V ^ ds ^J 



t > 0. 



Rearrangement therefore yields 



poo ( n \ 

J Q ( l|V«||| + ^||VCi|||J(is<oo. (4.6) 



From (|3.37p there is a constant C > such that 



di" 

For e > there is T e > such that 

||Vu|||ds < e 

by IjOj) . With /(t) := ||Vu|||, s(t) := C||Vw|||, fc(t) := C, oi = e, a 2 = Ce 
and 03 = Cr application of Lemma I3.5I yields 

||Vu(t + r)\\l < (e/r + Cr) exp(Ce), t > T £ . 
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So choosing first r > and then e > small enough shows ||Vu(t)||| < 5 
for t > T £ + r =: T$. Note that from Remark 1101 the uniform bounds on (f> 
and V0, and (|3.54p we have exactly the same situation for V£i as for Vu. 
So the assertion for VQ is obtained in the same way. □ 

Proof, (of Theorem \1.3\) . Existence. 

Let (u°,c°) G lf 1 ' 2 (fi)x(lf 1 ' 2 (H)n J L 00 (fi)) withe? > and f n c°(x)dx = 
rrii > 0, let (u, c) be the corresponding global strong solution from Theo- 
rem [L2] and let (ifc)fc be a sequence with tk — > oo. The estimates derived in 
Section [3] show that for the trajectory (u, c) there is a constant C such that 

||«(*fc)||2 + ||V«(* fc )|| 2 + \\c{t k )\\ 2 + ||VCi(t fc )|| 2 < C 

independently of time, cf. Lemmas 13.81 13. 9| 13.121 and 13.161 So (u(tk),c(tk)) 
is bounded in W '^(Q) x PF 1,2 (f2). By the compact embedding W 1,2 (U) ^^ 
L 2 (Q), the trajectory is relatively compact in L 2 (Q,), so there is a subse- 
quence of (ifc)fc<=N> which will be denoted in the same way, and an element 
(u 00 ^ 00 ) G L 2 (0) x L 2 (n) such that 

(«(t*),c(t fc )) -> KV 00 ) inL 2 . 

This convergence implies convergence almost everywhere in Q of a subse- 
quence, thus by nonnegativity of Cj(ifc) we deduce nonnegativity of cf 3 . Since 
||ci(£fc)||i = m i w e readily approve ||c?°||i = m;, i = 1, . . . , N. The conver- 
gence of Cj(tfc) in L 2 implies convergence of <fi(tk) —> 4>°° in W 2,2 iyi) where 
denotes the unique solution to 

N 

1=1 

With Poincare's inequality we obtain 

IKtfc)||2<c||Vu(^)|| 2 ->o, 

due to Lemma H31 so u°° = 0. 

Recall the definition Q = Ci exp(zi(f>) . From the uniform L°°-bounds on 
(j) and c, cf. Lemmas 13. 101 and 13.131 we see that the sequence (Ci is 
bounded in L 2 (Q). Since VCi(tfc) — > in L 2 (0), cf. Lemma hA\ we infer 
that Ci is bounded in W^ 2 (Vt). Since the embedding W^ 2 {9) ^ L 2 ($7) is 
compact there exists C°° £ L 2 (0) such that, up to a subsequence, £i(*fc) — ^ 
C°° in L 2 (0). Because VQ(t k ) -> we deduce that Cf° 

is in fact constant 

and that Ci(ife) -> Cf° in W 1>2 (n). 

From the uniform L°°-bounds of <p(tk) and Cj(ifc) it follows that 
X)j ll c j°lloo + 1 1 0°° I |oo < C. Thus from the continuity of the function 
(x, y) i — y xexp(ziy) we have £?° = c°° exp(zj^ 00 ) almost everywhere. 

Hence if we can show that c°° £ W 1,2 (J1), then Lemmas 14.21 and 14.31 imply 
that (u 00 ^ 00 ) is a solution to (P s ). 
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Note that from the definition of Q we have 

exp(-Zj0)VCi = Vq + ZiCiVtp. (4.7) 

Since \7Q(t) — > in L 2 (U) as t — > oo and the sequence 0(t&) is uniformly 
bounded in L°°(Q), for the proof of Vc°° € L 2 it is sufficient to show the 
convergence of Cj(tfe)V <j){tk) in L 2 . We compute 

\\ci(t k )v4>{t k )-cfvn 

< \\a(t k ) - cfHaHV^IU + Hcf 3 HoollV^tfc) - V0°°|| 2 -> 0, 

as — )• oo, where we use the convergence of — > <^°° in W 1,2 (fi) and 
the convergence of Cj(ifc) — >■ c?° in L 2 (f2). So by (|4.7|) we see that Vcj(tfc) 
converges in L 2 (0) as well, therefore cf° G W 1,2 (fi) and (u 00 ,*^ 00 ) = (0,c°°) 
is a steady state to problem (P). Therefore c°° G VF 2,2 (r2) is a solution 
to (Pj). 

Note that the assertion concerning the pressure 7r°° is already contained 
in Lemma 14.31 
Uniqueness. 

Due to the fact that u°° = it is sufficient to prove that c°° with cor- 
responding (p°° is the unique nonnegative solution to (Pg). So suppose 
there are two nonnegative solutions c 00 ,^ 00 with f n c°°dx = f n cfdx = mi, 
i = 1, . . . , N. Note that j n cf 3 — cf°dx = mi — mi = 0. Let us denote 
the corresponding potentials by <fp° , (j>°° . From Lemma 14.31 we have Cf° = 
c?° exp(zj0°°) = const and (f° = cf 3 exp(zi<f>°°) = const. Since C^ 00 , Cf° > 0, 
we can take the logarithm. It holds log£?° = logc?° + Zi4>°° = const and 
correspondingly for Therefore, integrating by parts we obtain 

N 

= E / ( c " - £ D(log C - log cT + ^(0°° - <A°°))dx 
i=i ^ 



V 



^(tf - ar)W»°° - + - cT) log 4- dx 



i 

.oo Too\i|2 i HJ.OO loo II 2 



£ || V (A 00 -* 00 )\i + er 



+ / E( c ~-oiog^> , 



so c°° = c°° and therefore </>°° = (f>°°. 
Convergence 

It was shown that semi-orbits (u, c) are compact and any accumulation 
point (u 00 ^ 00 ) represents a steady state. The convergence (u(tk),c(tk)) — > 
(u 00 ^ 00 ) in this situation means convergence with respect to the topology 
in W ' a (Q) x W 1,2 (Q). From the existence and uniqueness of steady states 
and the fact that the sequence (tk)k£N was arbitrary, the convergence (|1.38|) 
follows. □ 
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5. Asymptotic behaviour 

Let U C R 2 be bounded and smooth, let (u°, c°) E W$(Q) x (W 1,2 (I2) n 
L°°(fi)) with c° > and corresponding masses mi = f n c®dx > 0, i = 
1, . . . , N , be given and let c°° G W 2 ' 2 (0) denote the unique solution to p^y) 
subject to j n cfdx = mj, i = 1, . . . ,N from Theorem 11.31 

Let us define 



S ■= Wffin) x lv G M/ 1,2 (0) n :Vi>0,J Vi(x)d, 

and introduce the function ^ : 5 — >• R, given by 



x = vrii 



^{u, c) :=E(u, c) - E{u°°, c°°) 

A? 



dMli + E / [ciloga-c™ log c°°)dx 
1=1 J Q 

+ f (IIV0H! - ||V^|| 2 ) + y (||^||I ian - H^lllao 

N r / „ \ 

oo\||2 



^HI + E X ( Q(log ^ - 1} + c *°°) dx + l l|v(0 

+ e / V^°° • V((p- (f)°°)dx 
Jn 



2 



/ - <£°°)dz + V / (ci - c°°)(l + log tfjdx. 

Using integration by parts and the fact that c?° exp(zj^°°) is constant we 
obtain 

e [ V^°°- V^-^dx + er [ (ffty-^dx 
Jn Jan 

N 

+ E / {ci-cT){l + \ogc?)dx 

i=l ^ 

JV 
N 

= / E^- c ")( 1 + 1 °g^ oo )^ = o. 



i=l 

Therefore 



*(«, c) = i|H| 2 + £ jf (c t (log ^ - 1) + c°°) 



(5.1) 



+ + Lilian 
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Let (u,c) be the global strong solution to (P) from Theorem 11.21 with initial 
data (it°,c°). Note that due to Theorem 1 1 . 1 1 and Corollary [TJ for almost all 
t > we have (u(t),c(t)) G S, hence t i-> ^(u(t), c(t)) is well-defined. The 
value $f(u(t), c(t)) represents the difference between the energy at some time 
i and the energy in the equilibrium state. From the monotone behaviour 
of V along solutions of (P), cf. Lemma 13.71 we see that \& is monotonely 
decreasing as well with the same dissipation rate given in (|4.4p . 

Proposition 1. Let a, b > 0. 

(a) It holds 

(y/a - Vb) 2 < a(\og a - log b) + b - a. 
This inequality remains valid in the limits a — > and b — > as well. 

(b) Letb>b >0 and K := l/b . Then 

a (log a - log b) + b - a < K(a - b) 2 . 

Proposition[TJcan be proven by elementary calculus employing the relation 

2 log a — log 6 1 . . 

< — r 2-<-=, a>6>0. (5.2) 



a + b a-b Jab 



Theorem 5.1. Let VL C M 2 be bounded and smooth, let (tt°,c°) £ W ' a (£l) x 
(W 1,2 (0) n L°°(Q)) with initial masses mi = f Q c°dx > 0, i = 1,...,N, 
be given, let (u,c) be the global strong solution to (P) from Theorem 
with initial data (u°,c°), and let \& be defined as in 115.1]) . Then there are 
constants C, uj > depending on the initial data such that 

¥(u(t),c(t)) < Ce- Wt for a.a. t > 0. 

Proof. Recall the definition of V in ([331) and note that ft V = Hi*^ - ): c (-)), 
so from (|4.4p we have 

^(u(t),c(t)) = -D(u(t),c(t)), (5.3) 

where 

* /■ D 

D(u,c) = \\Vu(t)\\ 2 + Y] / ^exp(-2^))|V0W| 2 dx 
^•/n Ci{t) 

with £j = Cjexp(zji^). The idea of the proof is to estimate ^ suitably in 
terms of the dissipation rate D and obtain exponential decay via Gronwall. 

Let us define the auxiliary variable 



*(«)== i Miexp^M^L]-!. (5.4) 
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Recall that c°°exp(z^ 00 ) = (°° > due to Lemma [Ol For a / 0, the 
gradient of tpi satisfies 

|V^(t)| 2 = ^77^ T77^ 1VCi(t) + Wt)V<Ki)| 2 

4 Ci (t)c°° exp(zi(0°° - 0(i))) 

rD (5-5) 
<— ^exp(-^(t))|V0(t)| 2 , 

where we use the fact that c?° exp(zj<^°°) = const and ||^||oo < C uniformly 
in time, so 

l|V«(t)||l + W^Mml < CD(u(t),c(t)). 
i=i 

For proving the assertion it is sufficient to estimate 

*(u(t),c(t)) < C7 ^||V«(t)|| 2 + ^||VVi(t)||ij (5.6) 

with some constant C > independent of t, since then, by (|5.3|) and Gron- 
wall, we obtain exponential decay of ^f(u(-), c(-)). Note that we only require 
(|5.6p to hold true for those states (u(t),c(t)) which are attained by the so- 
lution (it, c). 

By Poincare's inequality we already know that 1 1 1* 1 1 2 < C||Vit||2 for all 
u G W Q ' a (Q). So for (|5.6p it suffices to show that for given R,R'>0 we 
have 

*(c) := 

:=£/ (cQog 
i=l ^ V 

<cj^ 11 Willi 

(5.7) 

for all c G Z^r/ with 

z^jy := jue ^' 2 (fi)nr(si) : 

Vi>0,J Vi = mi, iff (v) < R, \\v\loo < R'\, 
where G VF 3 ' 2 (0) is the unique solution to 

N 

—eAcj) = Z{Ci in f2, 9„</> + t</> = £ on c%7. 
i=i 

Note that for given initial data with initial masses m; > the phase space 
of c is contained in Zrri for sufficiently large R, R' > 0, cf. Corollary [IJ 



AT 

1) + cf\ dx + |||V(0 - + y ||0 - ^11^ 
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Remark [U and Lemma 13.131 From the Sobolev embedding W 2,2 (Q) 
BUC(Q) we compute 

IMloo < C\\<t>\\w2,2 < C'\\c\\ 2 < C'Hclloo < C"B! (5.8) 

for c £ Zrr', so we observe that the corresponding potentials </> are uni- 
formly bounded in L°°(f2). 

For the proof of relation (|5.7p observe that by the boundedness of c?° from 
below due to Lemma 14.31 and Proposition Q] there are constants C\ , C 2 > 
such that 

Ci ( E II - v^Flli + 11^ - <miW) < *( c ) 

i 

^cifeiict-criii+ii^-^ii^) 

i 

(5.9) 

for all c > and corresponding potentials <j>. 

We will derive (|5.7p indirectly. So suppose relation (j5.7j) does not hold. 
Step 1. 

Then there is a sequence — y 00 and (c fc )fc C ^,_r' with corresponding 
solution to 

-eA(f> k = E Zi c i\ m ^ <^ fc + T <^ = C, ° n dft, 
and auxiliary functions ip k , given by 

such that 

l k J2\\Wi\\ 2 2<y(c k )<R- (5-11) 

i 

Note that for c k , ^ superscripts refer to the index of the sequence and 
subscripts to the corresponding components. Without loss of generality we 
may exclude the case that c k = c°° for some k £ N, since in this case 
inequality (j5.7|) is trivially satisfied. 

By relation (|5.11|) we see Vip k — > in L 2 (Q.) for k — > 00. Due to the 
boundeds for c k and (j) k , implied by the definiton of Zr^ri and (|5.8p . we 
see that \\tjj k \\ 2 can be bounded by a constant independent of k, cf. ()5.10p . 
Hence the sequence (tpf)k is bounded in W l,2 (Q), and therefore ^ — 1 ipi in 
VF 1,2 (f2) and ^* — > ipi in L 2 (f2) via compactness. Because both i/)f — > ipi 
and Vt/jf — > in L 2 (S7), it follows that 

ip k -> $j in W 1,2 («) as fc -»• 00, (5.12) 

where VV>« = 0, hence ^ = const. 
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By dSSJ we have - <f>°°\\w^ < C, so there is 4> G W x ' 2 {n) such that 
cj) k -± $ in W 1 ' 2 ^), hence c/> fc -> 4> in L 2 (ft), as fc -> oo, (5.13) 
again using compactness. From the uniform boundedness of ||</> fc ||oo by (|5.8p 



and its convergence almost everywhere we deduce ||<^||oo < oo. Let 

6i := c°°(l + V^) 2 exp(^(0°° - 0)). (5.14) 

By the definition of ^ in f|5. 10|) we have = c°°(l + V'f) 2 exp(zj(^ 00 — (j) k )). 
Straight forward manipulations reveal 

c k -di = c? exp(*(0°° - <p k ))[(l + ^ fe ) 2 - (1 + ^) 2 exp(%(0* - 0))] 
= ef> exp(*(0°° - fe ))[(^ " ^) 2 + 2(Vf - + m4 ~ i>i) 
+ (1 + ^) 2 (l-exp^ -<£)))]. 

Using the a priori estimates ||</> fc ||oo + Halloo + Halloo + ||c°°||oo < C and the 
fact that VF 1 ' 2 (0) «^ L 4 (Q) gives 

||cf-Q|| 2 <C(||^-^|| 2 + ||^-^|| 2 + ||^-<A|| 2 ) ^Oas fe^oo. (5.15) 

By continuity it is true that — eA(j) = Zi&i in £1, and d u (p + r<j) = £ on dfi, 
and 

Cj exp(zi4>) = c°° exp(zi<p co )(l + V^) 2 = const. 

From relation (|5.14p we readily approve that Cj G VF 1,2 (il), since both ipi 
and c?° exp(zj(/) 00 ) are constant and 4> £ P^ 1,2 (0). Thus (c,<f>) is the unique 
steady state to (P) from Lemma [4.3l and Theorem ll.3| i.e. (c, 0) = (c°°, </>°°). 
So relation ()5. 15[) shows 

c k -»• c°° in L 2 (0) as A: ^ oo. (5.16) 

For the sequence (^(c fc ))fc this implies 

Afe := *(c fc ) -> (5.17) 

by 

Step 2. 

Recall that we assume / c°° for A; G N, so A 2 , = > for all A; G N. 

Having obtained the convergence of the sequences {ip k )k, {4 )k )k-, (c k )k and 
Afc in (|5.12p . (|5.13p . ([5.16P and (|5.17p we define new variables 



«?:=f, rf^ 5 ^, X fc :=^^. (5.18) 

Again for uj^ and y k superscripts represent the index of the sequence and 
subscripts denote the respective components. Note that for k G N the 
function x k is the unique solution to 

- eA X k = z iVi in n . and 9 ^ + rx* = on dSl. (5.19) 



A- 



t°°||2 



58 D. BOTHE, A. FISCHER AND J. SAAL 

Dividing (|5TT|) by l k \\ it follows that ]T\ ||V^ fc ||| < l/l k , so Vw k -> 
in L 2 (£l). We shall argue in an analogue way as for (|5,12p . For the 1?- 
boundedness of w k , observe that ip£ can be rewritten as 



With the help of (|5,9p we conclude 

N , N 

E n^fni + ii** - <t>°°Ww^ < c( E ii - v^f Hi + 

i=l ^ i=l 

< C'^{c k ). 

(5.20) 

Dividing (|5.20p by A? = ^(c k ) we have \\wf W2 < C" independent of k, hence 
w k m in W l ' 2 (Q) as 00, (5.21) 

where u) is constant. Also from the division of (|5.20p by A| we infer in the 
same way as for (|5.13p the existence of a \ G W l ' 2 {£l) such that 

X k x m W 1,2 (fi), hence x fe -»• X in L 2 (fi), as A; -»• 00. (5.22) 
With relations (lOTj) and (px22l) we define 

&:=C? (2«; i -^x)- (5-23) 
For the difference y\ — yi we obtain 



X - ((1 + exp(^°° - cf> k )) - l) - 2tB 4 + za) 
f (2|i(exp(^(0- - - 1) +2 (|L _ ^ + Zi U _ 



A fc Afc 
= C °°(2^(exp(^(^ - fc )) - 1) + 2K fc - t&O + z t (x ~ X k ) 

+ X k (w k ) 2 exp( Zl (^ - fe )) + ^^( X fe ) 2 exp(^(^ - fe ))) 
for some function h € [0, 1] using Taylor's formula. Hence 

U - mh < C(H k - <P°°\\l + \\wf - mh + \ k + \\x h - zh) ->• 0. 

By continuity it follows from (|5.19p that 

— eAx = E z iVi ™ ^' an( i ^vX + T X = on <9$7, (5-24) 
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and that \ k ~ t X m W 2,2 (£l). Therefore, using integration by parts and 
representation (I5.23p . we obtain 

0<e||Vx||l + er||x||i |an + / £ <?(2wi - z t xf 

Jn i 

= -e [ xA X + I y j cT{2w l -z l x? 
Jn Jn i 

= T2 z iXlli + m{2wi - ZiX) = 2 V] / Vim = 0, 

; Jtt i J& 

since JnVi = for all k. Therefore x = 0, u) = and finally = 0. 
However, dividing (|5.9p by A? now implies that 

1 < Yl \\M + llx*llv^ -»• as fc -> oo, 

j 

a contradiction. □ 

Now the proof of Theorem 11.41 remains an easy task. 

Proof, (of Theorem \1.4\ l- From the definition of *5f in f)5. 1[) it is clear that 
^ > 0, so Theorem 15.11 implies ||u(i)||2 < Ce""*. The claim for c follows 
with the help of the estimates in Lemma 13.131 Theorem 15.11 and relation 
(lEl by 

\\<H -cr\\x < C||(V* + V^)(v^- V^)ll2 < C'^Cu.c) 1 ^ < c'^-"'*. 

□ 
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